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SUMMARY 


This is the final report on NASA-MSFC Contract No. 
NAS8- 35594. The objective of this work was to improve our 
understanding of the variation of the aerosol backscattering at 
10.6 within the free troposphere and to develop a model to 
describe this. The analysis combines theoretical modeling with 
the results contained within three independent data sets. The 
data sets used are those obtained by the SAGE I/SAM II satellite 
experiments, the GAMETAG flight series and by direct backscatter 
measurements. The theoretical work includes use of a bimodal, 
two component aerosol model, and the study of the microphysical 
and associated optical changes occuring within an aerosol plume. 
A consistent picture is obtained, which describes the variation of 
the aerosol backscattering function in the free troposphere with 
altitude, latitude, and season. Most data are available and 
greatest consistency is found inside the northern hemisphere. 
The scarcity of data and lesser agreement within the southern 
hemisphere indicates the need for further experimental and 
theoretical studies. 
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1. INTRODUCTION 


The increased interest, in recent years, in the development 

and deployment of a Global Wind Measurement Satellite System 

(WINDSAT) (NOAA 1981), has led to an awareness of our lack of 

knowledge of the global characteristics of the free tropospheric 

aerosol. WINDSAT, as presently proposed, would use CO^ doppler 

shift lidar to measure wind velocity from the boundary layer to 

the lower stratosphere. It would rely on backscattering from the 

atmospehric aerosol for its signal and thus a knowledge of the 

global behavior of the aerosol backscattering function, 8 , is 

co 2 

essential for the system design. In a previous study (Deepak et 

al., 1982), a broad survey was made of the problems associated 

with modeling the aerosol backscattering function and the 

principal characteristics of its global variation were 

determined. These included a decrease in fi from the boundary 

co 2 

layer to the upper free troposphere and a strong latitudinal 
gradient with minimum values occurring over the southern oceans. 


The present study focusses on the use of three data sets to 


develop a global model for g 


CO^ 


, including variation with 


latitude, altitude, and season. Although present experimental 
work is exploring the use of alternative wavelengths between 9 
and 11 pm, the work in this report concentrates on the use of the 
standard wavelength, X = 10.6 pm. The data sets will include 
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those obtained from the Stratospheric Aerosol and Gas Experiment 
I (SAGE I) and Stratospheric Aerosol Measurement II (SAM II) 
satellite systems, the GAMETAG experimental series, and C0 2 lidar 
measurements. Individually, each data set has important 
limitations as far as the present objectives are concerned. 
Viewed collectively, as a data set which must be made 
self-consistent, they enable a significant advance to be made on 
previous global models. The SAGE I/SAM II data set is the result 
of a direct global measurement of aerosol extinction at a 
wavelength of 1 pm, a 1 0Q (McCormick et al., 1979). The 
satellites were designed for measurement of the stratospheric 
aerosol extinction but significant penetration of the troposphere 
down to an altitude of about 5 km is found, enabling the 
latitudinal and seasonal characteristics of the upper free 


tropospheric aerosol 

to 

be mapped. 

This 

information 

is for a 

wavelength of 1 pm 

and 

to obtain 

the 

equivalent 

map for the 

aerosol backscatter 

function at 

10.6 

e 10 . 6 . 

assumptions 


concerning the aerosol size distribution and composition must be 
made so that a conversion factor can be calculated. The GAMETAG 
data set consists of previously unanalyzed data obtained during 
two latitudinal flight surveys over the Pacific Ocean in 1977 and 
1978 (Patterson et al., 1980). Its analysis has been carried out 
under subcontract by the Georgia Institute of Technology under 
the direction of Dr. E. Patterson. These surveys measured the 
aerosol size distribution and composition and the data has been 
used to calculate the aerosol optical properties at 10.6 pm and 


i 
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other wavelengths. Its limitations are in terms of the small 
range in flight altitude (boundary layer or 5-7 km for the 
majority of the data) and the fact that these flights represent 
single time and space cross-sections. Uncertainties also exist 
because of the limited sampling volumes of the instruments used, 
resulting in significant statistical fluctuations in some of the 
particle counts. The direct lidar measurements of ^ appear to 
offer the most accurate input into a global model. 
Unfortunately, the amount of published data is very small; this 
situation may be expected to improve as several groups are 
currently making or planning to make measurements. All of the 
published data is confined to the northern hemisphere, and, 
although experiments are planned for the southern hemisphere, 
data for that hemisphere are not yet available. In addition, all 
these measurements have been made during or following a period of 
significant volcanic activity. Since 1979, several volcanoes 
have injected material into the stratosphere (Kent & McCormick, 
1984) and the recent eruption of El Chichon (McCormick et al., 
1984) is still profoundly influencing the stratospheric aerosol 
content. As will be shown later in this report, this volcanic 
influence extends into the upper troposphere as well as the 
stratosphere. WINDSAT cannot rely upon volcanic aerosol for its 
operation and present CC> 2 lidar measurements are thus not 
representative of the conditions under which it may have to 
operate. 
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In addition to the use of the above data sets, this report 
contains the results of further theoretical studies. In the 
first of these, modeling has been carried out of the 
microphysical processes and changes in optical properties 
occurring in long-distance transport of an aerosol plume. The 
aerosol in the free troposphere over the remote oceans originates 
mainly from long-distance transport of aerosol injected over the 
continents. Changes in the aerosol size distribution, 
particularly for the larger particles which sediment, occur 
during this transport, resulting in corresponding changes in 
g 1 . The other modeling which has been done is of an improved 
optical model in which a bimodal aerosol size distribution is 
used in conjunction with a two component aerosol. The main 
application of this model has been to the calculation of 
conversion factors to derive 8 , „ ^ from the SAGE I/SAM II 
extinction at 1 pm. 

The organization of this report is as follows: Sections 2 
and 3 present, respectively, the results of the analysis of the 
SAGE I/SAM II and GAMETAG data sets. As the latter contains a 
large number of plots representing the results of the analysis 
for individual stages of the 1977/19-78 flights, these have been 
removed from Section 3 and placed in the Appendix to this 
report. Section 4 contains the results of the microphysical 
modeling study and Section 5 presents the results of the bimodal 
size distribution modeling, and its use in the conversion of the 
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SAGE/SAM XI 1 pm extinction to the 10.6 pm backscattering 
function. Section .6 describes the available C0 2 lidar data set 
and Section 7 is an intercomparison of the results of the 
previous sections. Section 8 summarizes the main results, 
indicating the areas where further development is needed. 


J 
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2. ANALYSIS OF THE SAGE I/SAM II DATA SET 
2.1 SAGE I AND SAM II TROPOSPHERIC OBSERVATIONS 

The SAGE I (Stratospheric Aerosol and Gas Experiment I) and 
SAM II (Stratospheric Aerosol Measurement II) satellite 
experiments contain sun photometers designed to measure the 
extinction produced by stratospheric aerosols at a wavelength of 
1 pm (McCormick et al., 1979). SAM II r which was launched in 
October 1978, and which is still operational, consists of a 
single channel sun photometer only. The orbit of SAM II is such 
that its observations, which are made at satellite sunrise and 
sunset, occur between latitudes of 64 S and 81 S and 65 N and 83 
N, two measurements being made on each orbit. SAGE I commenced 
its observations in February 1979 and, owing to a faulty 
satellite power supply, these were terminated in November 1981. 
In contrast to SAM II, the SAGE I coverage is nearly global, the 
latitude of observation moving during a six-week period through 
about 120° of latitude, the latitude extremes being approximately 
74° N & S. SAGE I also contains an aerosol channel at a 
wavelength of 0.45 pm in addition to the 1 pm channel. Useful 
data on the former channel is limited to altitudes above 10 km 
and, in this report, our analysis has been confined to the 1 pm 
data only. The SAGE I data set is also limited by the fact that 
observations were confined to sunsets only, after the first few 
months of observation, in order to conserve satellite power. 

As noted above, SAGE I and SAM II were designed for the 
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measurement of stratospheric aerosols and it was anticipated that 
tropospheric measurements would be hindered or possibly prevented 
by the presence of high altitude cloud. The former is indeed the 
case but significant tropospheric penetration does occur. An 
example of a SAGE I profile showing penetration down to an 
altitude of 6 km is shown in Fig. 2.1. (a). It may be noted that 
the aerosol extinction increases smoothly with decreasing 
altitude, no apparent discontinuity occurring at the tropopause. 
Fig. 2.1(b) shows a different profile with a strong enhancement 
at an altitude of 10 km, presumed to be due to high altitude 
cloud. Table 2.1 shows the total number of observations made by 
SAGE I over a 3-month period and the frequency of penetration to 
various tropospheric altitudes for different latitude bands. It 
can be seen that good penetration (> 50%) is obtained down to an 
altitude of 8 kilometers or less. Greatest frequency of 
penetration occurs in the southern hemisphere and, as might be 
expected, least occurs in the equatorial zone with its higher 
tropopause. This table shows that one might expect to obtain a 
resonable description of the free troposphere aerosol, on an 
average basis, down to an altitude of perhaps 5 or 6 km. Below 
this altitude, individual profiles, when available, will be 
accurate but average values, which represent only 20% - 30% of 
the total potential data set, may be expected to be biased toward 
more transparent atmospheres and lower extinction. 



Table 2.1. Relative Frequency of SAGE Observations in the 
Troposphere (March - May 1979) 


Altitude 

(km) 


Relative Frequency (%) 


60°S - 20°S 

20°S - 20°N 

20°N - 60°N 

20 

100 

100 

100 

18 

100 

96 

100 

16 

100 

78 

99 

14 

98 

61 

98 

12 

89 

53 

89 

10 

73 

48 

68 

8 

58 

45 

52 

6 

44 

36 

34 

4 

30 

17 

15 

2 

6 

0 

4 

TOTAL 

OBSERVATIONS 827 

290 

455 





CM 



rM 



(i/w) ganimv 
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(a) No high-altitude cloud present. (b) High-altitude cloud present. 

Figure 2.1. SAGE 1 jjm extinction profiles showing penetration into the troposphere. 





2.2 EXTINCTION PROBABILITY DISTRIBUTIONS AND USE OF THE MEDIAN 
EXTINCTION 

Examination of the SAGE I and SAM II data shows that 
whereas, in the stratosphere, extinction values at a given 
altitude and latitude are fairly tightly concentrated about a 
mean level, the same is not true in the troposphere. In the 
latter case, extinction values may vary over several orders of 
magnitude, the higher values probably being due to attenuation by 
thin (sub-visible) cloud. Under such conditions, the use of the 
mean extinction as a measure of central tendency is of doubtful 
value. In order to obtain a better description, we have 
examined, in some detail, the extinction probability 

distribution. 

Figure 2.2(a) shows the result of binning three months of 
SAGE I data by altitude and extinction. The data has been taken 
over the latitude band 20° - 60° N, between March and May 1979, 

and the numbers on the diagram represent the numbers of 

- 5 -1 

observations falling into each (1 km) x (2 x 10 km ), 

altitude-extinction bin. In the stratosphere, the extinction 

values are tightly bunched at a given altitude and there is no 

problem about defining and using a mean level. In contrast, in 

the troposphere below about 12 km, the values are widely 

separated and there is a significant number of very high values 

-3 -1 

(Extinction >10 km ) which are shown in the bins at the 

bottom of the diagram. These values probably represent 
attenuation by sub-visible cloud (Rao, 1975; Uthe & Russell, 
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1977). Superimposed on^ t?he diagram are lines showing the 
cumulative probability levels for the extinction. In the 
stratosphere they are close together, in the troposphere, they 
are quite widely separated. Rather than averaging the extinction 
values at a given altitude to form a mean value, we have chosen 
to use the 50% probability level, or the median, as our measure 
of central tendency. In the stratosphere, it approximates very 
well to the mean level. In the troposphere, it defines an 
aerosol level which is both useful and meaningful. The numerical 
value of the median is not appreciably affected by the inclusion 
of very high extinctions due to cloud and its relationship to the 
extinction probability distribution makes it compatible with 
published CO lidar backscatter observations, e.g. (Post et al., 
1982). A second example of binnned aerosol data and the 
associated probability levels is shown in Fig. 2.2(b). This is 
for the same time period as Fig. 2.2(a) and for the same latitude 
band in the southern hemisphere. The marked difference between 
the two diagrams with less extinction occurring in the southern 
hemisphere is a systematic feature of the SAGE I and SAM II 
data . 

Figure 2.3(a-e) shows cumulative extinction probability 
distributions for five latitude bins and various altitudes. The 
data is for the period March-May 1979 and is for both SAGE I and 
SAM II. The choice of axes and scales on these figures is such 
that log normal probability distributions appear as straight 
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Figure 2.2. Scatter plots showing the distribution of SAGE 1 fjm extinction values at different 
altitudes. Superimposed on the scatter plot are lines showing the 20%, 50%, and 
80% cumulative probability levels for the extinction. 



Altitude (km) 


,1 


« p ORIGINAL PAGE IS 
OF poor QUALrrr 



o 


« 

*•» 


* 

o 


« « «« rTriMi f*cr oc • 9eo 

I i I i i i » i « i i • • i * i i i i i i i i i i » i » i i i » « i i i » i » i i J * * i » J * • ' 

O ■* « « ON#9COA.-fOCON<f«ftCN<«KQN9«aON<t«<OM9 dCCN4««CN4ei 


o o o 

— CM 


o o o 

IO 9 IO 


o o o o o 

<0 h- CO 0> O 


( i-UJ^t g_OI) U0!PUUX3 


<u 

s . : 


cn 

hs 

05 


>« 

m 


x 

CJ 

s_ 

CO 


ro O) 
cn C 

11 
5 x 
> in 

C w 
o <U 
"43 .£ 

u — 

.£ a) 

•f !_ 

X (0 
0) 


0) 

E ° ~ 

r i° 


LU 
O co 


< 

to 


M- ® 

O 4= fc* 


(/) 


© 

co 


o 

CM 


<u 

a) 

x 

a 

cn 

E 

a) 

x 

c 

s_ 

D 

X 

■M 

3 

O 

(/) 


c •; 


4-> 


CO 

3 

T5 

X 

X 

a) 

o 


in 

s_ 

"si 

+■> 

cn 

O 

a 

a 


£ 

9) 

> 

<u 

X 

C 

a) 

a 

"43 

_co 

+J 

3 

3 

04 

c 

tn 

E 

3 

(J 

i 

in 

o\° 

o 

9) 

o 

X 

■o 

00 

cn 

3 


cn 

+■» 

+- T3 
"43 C 

o 

CO 

CD 

a 


-S 


4-* 

o\o 

s_ 

£ 

o 

(1) 

■(-< 

L. 

LD 

4J 

0 


co 

o 

4- 

It 

o\o 

O 

U) 


CN 

• 



CN 



CN 



a 



L. 



3 



04 


LL 




4 


2-8 



ALTITUDE (km) 


* * A 

' ; :::••>*< 


ID O T CD 


<i> $ <> * 


\ • 


% <> 

V o 

9 

\ o 
% $ 

♦J 


90 * 

o g 


(,-UJ5|) N0I10NI1X3 


« + 

■° 0 0 0 


0. s 0 

6 

0 0 * « 

0 9 

oO- 0 

■( 

o 

o 

o 

o 


v CPiv 

A « 

9 .<> 0 

E 

0 o 

UJ 

§ CD O ^ 00 
r 

0<> 

V* 

b 

99 « 

^ 4> «> <> ^ 



(i-W>0 N0I10NIIX3 


Figure 2.3. Cumulative probability distributions for SAGE I and SAM II extinction, 
March - May, 1979. 
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Figure 2.3. Cumulative probability distributions for SAGE 1 and SAM il extinction, 
March - May, 1979. 





lines. Although no detailed study has been made of these 
distributions, several interesting points emerge. 

(1) In the stratosphere, the probability distribution 
appears to be log-normal. The only exception to this is the 
distribution for an altitude of 18 km shown in Fig. 2.3(e). At 
this time of year, the stratospheric aerosol at this altitude is 
strongly differentiated depending upon whether it is observed 
inside or outside the polar vortex (McCormick et al., 1983). The 
two slope nature of the probability distribution reflects this 
division in the data, which includes measurements made both 
inside and outside the vortex. 

(2) Near the tropopause [e.g. , at altitudes 14 and 18 km in 
Fig. 2.3(c)], the distribution shows a strong discontinuity in 
slope. The lower section of the distribution is believed to be 
due to aerosol , the higher and steeper section to the presence of 
sub-visible cloud. 

(3) The steeper slope and vague linearity of the probability 
distributions at lower altitude in the troposphere indicates the 
wide range of aerosol conditions included in the data set. 

2.3 VARIATION WITH LATITUDE, ALTITUDE, AND SEASON 

In order to study the variation of the aerosol extinction 
with latitude and season, the SAGE I data was grouped into seven 
latitude bands and the SAM II data in four latitude bands, as 



shown in Table 2.2. Figure 2.4(a) shows the altitude variation of 
the median SAGE I extinction for the six latitude bands covered 
during the period March-May 1979. Apart from data in the 60 N - 
75 N latitude band, there is a general decrease in aerosol 
extinction with increasing altitude. In examining these and 
other data presented in this section, it should be remembered 
that below an altitude of 5 or 6 km, the fractional penetration 
is less than 50% and the data may have a systematic bias. This 
is particularly so when the extinction is high, as in the 60 N - 
75 N latitude band, and it is doubtful if in this case the 
decrease in median extinction for altitudes below 5 km is 
representative. It is more likely that observations are being 
made down to these altitudes only when the atmosphere is 
relatively clean of both aerosol and cloud. A secondary feature 
of the variation with altitude is the greater extinction observed 
in the upper free troposphere within the equatorial belt (20 S - 
20 N) as compared to the other latitude bands. This reflects the 
higher tropopause level (~16 km against ~ 12 km for mid-latitude) 
and must indicate the effects of convection in raising the 
aerosol to the higher levels. 

Apart from the variation with altitude, the most important 
feature is the marked latitude asymmetry. At an altitude of 6 
km, the variation in extinction between 60 N - 75 N and 40 - 60 S 
is approximately one order of magnitude. Although this asymmetry 
is present at all times of the year, it is at a maximum in 



TABLE 2.2. 


(a) SAGE I 


(b) SAM 11 


Latitude Bands Used in the Analysis of 
SAGE 1 and SAM 11 Data 


60 N - 75 N 
40 N - 60 N 
20 N - 40 N 
20 S - 20 N 
40 S - 20 S 
60 S - 40 S 
75 S - 60 S 


75 N - 90 N 
60 N - 75 N 
75 S - 60 S 
90 S - 75 S 
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Figure 2.4. Median 1 pm aerosol extinction profiles for SAGE I data shown as a function of 
latitude band. 



March-May, both in 1979 and in the following years. It seems 
possible that the occurrence of particularly strong extinction in 
high northern latitudes at this time of year may be connected 
with the observation of arctic haze (Schnell, 1984) at somewhat 
lower altitudes. 

Figures 2. 4(b) -(d) show the equivalent aerosol extinction 
plots for June-August 1979, September-November 1979, and December 
1979-February 1980. The features noted above are present in all 
these plots, the latitude asymmetry is, however, less in Figs. 
2.4(c) and (d). The peak in aerosol extinction visible in the 20 
S - 20 N latitude band in Fig. 2.4(d) at an altitude of 19 km is 
caused by the injection of material from the eruption of the 
Sierra Negra volcano on November 13 , 1979 (Kent & McCormick, 
1984). Figures 2.5(a) and (b) show the high latitude aerosol 
extinction as measured by SAM II during March-May 1979 and 
September-November 1979. The tropopause altitude is lower (~ 10 
km) and colder temperatures in the antarctic stratosphere have 
produced an enhancement in the median aerosol extinction. In the 
troposphere, the hemispheric asymmetry observed at lower 
latitudes is still very evident. 

Figures 2.4 and 2.5 may be used to study the seasonal 
changes which occur in the aerosol extinction. An alternative 
way of displaying this which shows the seasonal variation very 
clearly is given in Fig. 2.6(a) and (b). In these figures, the 
SAGE 1/SAM II extinction at an altitude of 6 km has been plotted 
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Figure 2.4. Median 1 jam aerosol extinction profiels for SAGE I data shown as a function of 
latitude band. 
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EXTINCTION 




against latitude for the same four seasons as shown in Fig. 2.4 
and 2.5. Both figures show the hemispheric asymmetry; in 
addition, they show a clear superimposed seasonal variation. In 
both hemispheres, maximum aerosol extinction is found in local 
Spring-Summer and minimum extinction in local Fall -Winter. As 
most aerosol at this altitude is derived from continental 
sources, the higher values in the northern hemisphere most likely 
reflect the greater land surface as compared to the southern 
hemisphere. The Spring-Summer maximum is most probably related 
to the increased convection over land at that time of year. 

Another way of displaying the data in Figs 2.4 and 2.5 is 
shown in Fig. 2.7(a) and (b). In these figures, the aerosol 
extinction is plotted in the form of a contour diagram as a 
function of latitude and altitude. These figures will later be 
converted to equivalent plots for the backscatter function at 
10.6 pm . 

2.4 VARIATION WITH SURFACE TYPE 

Most of the aerosol in the free troposphere over the remote 
oceans is derived, not by convection from the ocean surface, but 
by long-range transport of aerosol from over the continental land 
masses. Some modification of the aerosol optical properties 
might be expected to occur during this transport and in Section 4 
a microphysical model is presented which investigates this 
variation theoretically. In the present section, we report the 
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Figure 2.7. Contour plots of free tropospheric aerosol extinction at 1 pm from SAGE I and SAM il. 



results of an examination of the SAGE I data set for experimental 
evidence of such a variation. For the purposes of this study, 
the global surface has been divided into 10-degree 
latitude-longitude squares and each square categorized according 
to one of the five classes shown in Table 2.3. Two classes of 
remote ocean have been defined as it was felt that the antarctic 
continent was unlikely to be a major source of aerosols. A 
global map showing the classification of each 10-degree 
latitude-longitude square on the surface of the globe is given in 
Fig. 2.8. 


Each SAGE I observation has been classified according to the 
surface type beneath the observation position, as well as for 
season and latitude. Median extinction values have been 
calculated for each altitude. The analysis shows very little 
variation of aerosol extinction with sub-surface type, typical 
results of this analysis being shown in Fig. 2.9(a) and (b) . In 
these figures, it can be seen that any systematic difference 
between the extinction over land or ocean is less than the error 
in the median values as shown by the error bars in the figures. 
The only positive result to be obtained from the analysis is for 
the equatorial region (20 S - 20 N ). Figure 2.10(a) shows 
histograms of percentage departures from the zonal mean 
extinction at an altitude of 6 km for the four surface classes 
found in the equatorial region. Individual data points in the 
histograms are for the medians of 3-month data sets and the 



TABLE 2.3. Set of Surface Types Used 
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Figure 2.8. Division of the global surface area by surface type; 

L - Land R1 - Remote Ocean I 

C - Coastal R2 ~ Remote Ocean II 
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three-month median levels, the eleven 
points within each histogram covering 
the entire 33-month SAGE I data set. 



eleven points within each histogram represent the entire 33-month 
SAGE I operational period. The median values of all the data 
within each histogram shows a small but steady shift from +9.7% 
for land to -9.3% for remote ocean. This very small change is in 
the expected direction and is indeed comparable to the 
statistical errors. At other latitudes, no such systematic 
variation is observed and it is likely that this difference is 
related to the magnitude of the zonal wind velocity. Figure 
2.10(b) shows the time averaged zonal wind at the 500 mb level as 
a function of latitude (Lorenz, 1967). The lowest velocities are 
observed within the equatorial belt giving greater transport 
times for movement of aerosol from over land to over ocean. 
Corresponsively, greater changes in microphysical properties and 
related optical extinction are thus likely to occur. 

2.5 VOLCANIC EFFECTS 

It is well known that volcanic eruptions inject solid and 
gaseous materials into the stratosphere and that aerosols formed 
and deposited there have lifetime of months or years 
(Deirmendjian, 1973; Deepak, 1982; Kent and McCormick, 1984). 
Solid particles are also injected in the troposphere where the 
majority of them are presumed to be fairly quickly removed by 
sedimentation and washout. Examination of the SAGE I and SAM II 
data set show that the upper troposphere, as well as the 
stratosphere, has a long-lived enhancement in aerosol 
extinction. Examples of this enhancement are presented in Figs. 
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2.11 (a) thru (d), which shows vertical median extinction 
profiles for the four seasons between March 1980 and February 
1981. In 1979 the only volcanic eruption of significance was that 
of Sierra Negra on November 13, 1979, for which the major 
stratospheric effects were not felt until early in 1980. The 
stratosphere (and upper troposphere) in 1980 was affected not 
only by the after-effects of the Sierra Negra eruption but also 
by the eruptions of St. Helens (46 N, May 18, 1980) and Uluwun (5 
S, October 7, 1980). The 1980-81 seasonal extinction profiles 
(Fig. 2.11) may be compared with their equivalents for 1979 (Fig. 
2.4). Those for March thru May are similar, except for the 
stratospheric enhancement produced by Sierra Negra and visible in 
the 20 S - 20 N latitude belt in March thru May, 1980. Those for 
June thru August are very different. In 1980, the free 
troposphere between 40 N - 60 N and 60 N - 75 N is profoundly 
modified following the eruption of St. Helens. In the case of the 
60 N - 75 N band, the 1980 increase is visible down to an 
altitude of about 5 km, well beneath the tropopause altitude. 
Similar differences are observed within the same latitude bands 
for September-November and in the 40 - 60 N latitude band for 
December-February (no data is available for the 60 N - 75 N 
latitude band during northern winter); in all cases the 
enhancement appears to occur in the free troposphere down to an 
altitude of about 5 km. The stratospheric enhancement during 
this period and within these latitude bands is produced by the 
injection of material from the St. Helens eruption which was the 
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largest of the four eruptions discussed above. There seems 
little doubt, therefore, that the free tropospheric enhancement 
is related to the same source. It is noticeable that the 
enhancement reaches to maximum amplitude close to the tropopause 
( ~10 km at 60 N, see Fig. 2.7) and it is possible that the 
aerosol in the stratosphere is acting as a reservoir from which 
material is being fed into the upper troposphere by sedimentation 
and stratospheric-tropospheric exchange processes. The data 
shows many interesting and puzzling features which will be the 
object of future study. For example, the data shows that the 
maximum transfer of aerosol from the stratosphere to the 
troposphere occurs at high latitudes, whereas no such transfer is 
evident at low latitudes, although volcanic material is clearly 
present in the stratosphere (Fig. 2.11(a), (b) , and (d)). These 
differences may help resolve the relative importance of 
sedimentation, stratospheric-tropospheric exchange and general 
circulation as mechanisms for transfer of material from 
stratosphere to troposphere or vice-versa. 
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3. ANALYSIS OF THE GAMETAG DATA SET 


3.1 INTRODUCTION 

Section 3, along with Appendices A and B, are based on the 
Final Report submitted by Georgia Institute of Technology, 
Principal Investigator, Dr. E. M. Patterson, under IFAORS sub- 
contract F- 520 . 

The Global Atmospheric Measurement Experiment of 
Tropospheric Aerosols and Gases (GAMETAG) has as its objective 
the coordinated measurement of those atmospheric trace species 
(gases and aerosols) that are necessary to gain an understanding 
of atmospheric chemical properties and processes. In particular, 
the aerosol measurements of the phase I GAMETAG program were 
designed to measure the levels and types of tropospheric aerosols 
both in the planetary boundary layer and at mid-tropospheric 
altitudes in the free troposphere under a variety of remote 
conditions, to study aerosol chemical processes, and to study the 
optical effects of the aerosols as an indication of possible 
climatic impact. 

The aerosol measurements during the GAMETAG program have 
provided the most extensive data set available for 

mid-tropospheric remote area measurements. A detailed 
description of these aerosol measurements, emphasizing the 

aerosol chemistry, has been reported by Patterson et al. (1980), 
but only a limited analysis of the optical properties has been 
reported. 
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It is the purpose of this chapter to present in considerably 
greater detail some calculated visible and infrared optical 
properties of the measured aerosol in both the free troposphere 
and the boundary layer. The visible wavelength calculations were 
made for A. = 0.63 pm for comparison with on-board in-situ 
scattering instruments. Since one of the major aims of this work 
is the determination of aerosol properties in the 10 pm 
wavelength range for the analysis of the expected signal from an 
orbiting doppler lidar, the infrared calculations were made for A 
= 10.6 pm. Calculations were made of extinction and backscatter 
coefficients for each of these two wavelengths; ratios of 
backscatter to extinction were also calculated for each of these 
wavelengths . 

3.2 GAMETAG FLIGHT PROFILES 

There were two series of GAMETAG flights. One took place in 
August and September of 1977; the other occurred in April, May, 
and June of 1978. The GAMETAG flight patterns for the 1977 and 
1978 operations are shown in Figure 3.1. A listing of the data 
and location for each flight is given in Table 3.1; analysis was 
done for those days marked with +'s. Generally, the two flight 
series consisted of flights over continental and adjacent marine 
areas of North America and flights across the Pacific Ocean that 
were designed to cover as wide a latitude range as possible. In 
general, the flights consisted of sampling legs in both the 
mid-tropospheric region (5-6 km altitude) of the free troposphere 



6AMETAG FLIGHT OPERATIONS 



Figure 3.1. Flight tracks for the NCAR Electra aircraft 

during 1977 (solid line) and 1978 (dashed line). 
(From Patterson et al., 1980.) 
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TABLE 3.1. Detailed Flight Tracks During the 1977 and 
GAMETAG Flight Series 


DATE 



INITIATION 




1977 

Aug. 

7 

+ 

Denver 

Aug. 

8 

+ 

Portland 

Aug. 

9 


Anchorage 

Aug. 

11 


Anchorage 

Aug. 

12 


Churchill 

Aug. 

22 

+ 

Denver 

Aug. 

23 

+ 

San Francisco 

Aug. 

25 

+ 

Hilo, Hawaii 

Aug. 

26 

+ 

Johnston Island 

Aug. 

28 

+ 

Pago Pago 

Aug. 

31 

+ 

Pago Pago 

Sept. 

1 

+ 

Pago Pago 

Sept. 

2 

+ 

Johnston Island 

Sept. 

5 


Hilo, Hawaii 

Sept. 

6 


San Francisco 


1978 


April 

27 

+ 

Denver 

April 

28 


San Francisco 

May 

2 

+ 

Hilo, Hawaii 

May 

3 


Johnston Island 

May 

4 

+ 

Canton Island 

May 

6 


Fiji Islands 

May 

10 

+ 

Christchurch, N 

May 

11 

+ 

Christchurch, N 

May 

12 

+ 

Fiji Islands 

May 

13 


Canton Island 

May 

14 

+ 

Johnston Island 

May 

17 


Hilo, Hawaii 

May 

18 

+ 

San Francisco 

May 

27 

+ 

Denver 

May 

28 

+ 

Great Falls 

May 

30 


Whitehorse 

May 

31 

+ 

Whitehorse 

June 

1 

+ 

Great Falls 


1978 


DESTINATION 


Portland 
Anchorage 
Anchorage 
Churchill 
Denver 
San Francisco 
Hilo, Hawaii 
Johnston Island 
Pago Pago 
Pago Pago 
Pago Pago 
Johnston Island 
Hilo, Hawaii 
San Francisco 
Denver 


San Francisco 
Hilo, Hawaii 
Johnston Island 
Canton Island 
Fiji Island 
Christchurch, N.Z 
Christchurch, N.Z 
Fiji islands 
Canton Island 
Johnston Island 
Hilo, Hawaii 
San Francisco 
Denver 
Great Falls 
Whitehorse 
Whitehorse 
Great Falls 
Denver 
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(that portion of the troposphere not directly influenced by the 
surface) and in the planetary boundary layer (0-2 km altitude) 
with sounding data taken on the ascents and descents available to 
link the measurements in the two regions. The detailed flight 
tracks for each flight are shown in Appendix A. 

The data were taken during two seasons , with most of the 
flights centered around local noon for maximum photochemical 
activity. Because of the scientific objectives of the program, 
the flights were planned to avoid regions of strong convective 
activity and to avoid cloud penetration as much as possible. 
Thus the measurements emphasized clear air data. 

3.3 INSTRUMENTATION 

The size distributions discussed here were determined by 
means of single particle optical counters. As discussed in 
Patterson et al. (1980) the inferences from these optical 
particle counters were tested by comparison with filter and 
cascade impactor measurements . 

Two optical particle counters manufactured by Particle 
Measurements Systems (PMS) of Boulder, Colorado, were used on the 
National Center for Atmospheric Research (NCAR) Electra aircraft 
during GAMETAG. A combination active scattering aerosol 
spectrometer and classical scattering spectrometer probe 
(ASAS-CSSP) was mounted within the aircraft, and it sampled air 
from an aerosol inlet manifold, which also supplied the sample 
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air for the filter samples. The aerosol inlet manifold was 
designed to sample air isokinetically from the outside of the 
aircraft, to have a diffuser section to reduce the air speed 
within the sample tube, and to have a series of intakes for each 
of the aerosol measurement devices. All intakes were operated 
approximately isokinetically, and the entire system was 
wind-tunnel tested for sampling efficiency. Such testing showed 
sampling efficiencies of ~ 1 for particles with r 0.75pm. 
Sampling efficiency decreased for the larger particles, and so an 
externally mounted probe, the forward scatter spectrometer probe 
(FSSP), was used to measure the larger particles. Each of the 
probes has been extensively calibrated by the manufacturer and by 
the GAMETAG experimenters by using polystyrene latex spheres and 
glass spheres for the particle ranges used. Our determination of 
sizes from the optical particle counter data is based on the 
polystyrene and glass sphere calibration, with no explicit use of 
calculations for the differing response of the instruments to 
particles with differing refractive indices. Our calculations of 
the Mie scattering functions appropriate to the PMS optical 
geometry, as well as those reported by Pinnick and Auvermann 
(1979), indicate that for the range of sizes and refractive 
indices encountered in this experiment, the differing 
instrumental response characteristics will have no significant 
effect on our measured size distribution. 

The data from the PMS probes were transferred to the Electra 
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Electronic Data Management system tape every 2 s, but for 
analysis, these second-by-second averages were converted to 
longer time averages to reduce statistical fluctuations. One 
minute accumulations were used as the standard archived output 
from the particle probes. 

Although each probe is a multi range device (for a complete 
description see the ASAS and FSSP manuals, available from PMS), 
only one of the ranges for each probe was used in the analysis of 
the size data: the nominal 0.25 pm 4 : r 4 : 3.0 pm for the FSSP and 
nominal 0.15 pm < r 0.75 pm for the ASAS-CSSP. 

The response calculations for these optical particle 
counters exhibit double valuedness, and so the ASAS and the FSSP 
size data were combined into a smaller number of larger 
intervals. Nine combined size ranges were used in the 
calculations, six ASAS size ranges and three FSSP size ranges. 

The same ASAS size ranges were used for both the 1977 and 
the 1978 data sets. The FSSP' data for the two data sets was 
handled somewhat differently, however. For the 1977 data, the 
FSSP data was combined into three ranges consisting of Channel 1 
(FSSP Range 7), Channels 2-9 (FSSP Range 8 ), and the rest of 
the FSSP Channels (FSSP Range 9). Range 7, which overlapped the 
ASAS data, was not used and Range 7 consisted of Channels 2-5, 
Range 8 consisted of Channels 6-9, while Range 9 consisted of 
the rest. These ranges, size limits, and mean radii are shown in 
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Table I of Appendix B. 


3.4 DATA ANALYSIS 

The earlier analyses by Patterson et al showed that the 
typical size distribution exhibited a bimodal nature as shown in 
Fig. 3.2. Analysis of the filter and cascade impactor samples 
showed that these modes had different composition over land, with 
the submicron mode composed of sulfate and other secondary and 
combustion aerosols, and the supermicron mode composed of 
soil-derived aerosol particles. The marine free tropospheric 
aerosol was similar to the continental aerosol ; the marine 
boundary layer aerosol, by contrast, was dominated by a sea salt 
aerosol . 

Consequently, we have used more than one set of optical 
constants to describe the aerosol. Over continental areas, and 
in the marine free troposphere the particles with r > 0.5 pm were 
described by optical constants appropriate to soil aerosols, 
while the smaller particles with r < 0.5 pm were described by 
optical constants appropriate to ammonium sulfate. Optical 
constraints appropriate to a wetted sea salt aerosol were used 
for the total size distribution in the marine boundary layer. 
The optical constants used in these calculations are shown in 
Table 3.2. 

Calculations of the extinction coefficient a £ (expressed in 


3-8 



dN/dLOG ( R ) 


***Kno I I ertberg Doto*** 
DRTE=MRY 2. 1978 



STRRTING TIME =21:23: 0 
ENDING TIME =21:42: 0 


Radius (micron) 


Figure 3.2. A representative size distribution measured during 
the GAMETAG flights. 




Table 3.2. Optical Constants Used in Calculations 


Wavelength 

0.63 pm 

10. 

6 pm 

(NH 4 ) 2 S0 4 

1.5 - 0.005i 

1.98 - 

0 . 06i 

Soil-Aerosol 

1.9 - 0. Q05i 

1.74 - 

0.4i 

Sea Salt (Wetted) 

1.4 - O.Oi 

1.38 - 

0. 057i 
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-1 

terms of m ) and the backscatter coefficient 3 (in units of 
-1 -1 

m sr ) were made for 0.63 jam and 10.6 }im. A Mie scattering 
routine available at Georgia Tech was used for the calculations. 
The calculation procedure involved a calculation of the optical 
quantity of interest using appropriate optical constants for each 
size range. These average efficiency factors were multiplied by 
the average particle cross sectional area for the interval and by 
the number density of the particles to determine cr £ or 3 for each 
size interval. Values for the individual size intervals were 
then summed to determine a or b for the entire distribution. 

Calculations were made for the basic one-minute data set 

3 

that consists of counts/cm in each one of the nine channels. 
These one-minute averages were combined into five minute averages 
which are plotted, together with the standard deviation 
calculated for each set. These one-minute and five-minute data 
are plotted for each of the days analyzed. 

The time averaging corresponds to a spatial averaging along 
the aircraft flight .tracks as well. A one-minute average 
corresponds to a distance of ~ 9 km at high altitudes and ~ 7 km 
at low altitudes. Five minute averages correspond to a spatial 
averaging over ~ 50 km at high altitudes and ~35 km at low 

altitudes. 

In addition, longer time averaging was also done. The 
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one-minute particle averages were combined into 20 minute 
averages to produce a series of size distributions which are 
presented in Appendix B. The one-minute optical data were also 
averaged over flight legs or portions of flight legs (generally 
100-200 km) and plotted as a function of latitude to determine 
latitudinal profiles along our flight tracks. 

We will discuss the calculated optical properties in terms 
of the daily flight data and the latitude profiles. We will not 
further discuss the size distribution data. 

3.5 DAILY FLIGHT DATA 

Detailed optical calculations were made for twelve of the 
1978 flights and ten of the 1977 flights. These calculations 
were made as described above for the optical properties of 
interest: extinction coefficients at visible and infrared 
wavelengths, backscatter coefficients at these wavelengths, and 
backscatter-extinction ratios. These quantities have been 
plotted against time into the flight for each of the days 
analyzed. This data set is shown in Appendix A, arranged by 
flight days. In this Appendix, each data set consists of ten 
pages including a table and nine figures (Figures a-i) arranged 
as follows: 

1. An initial table giving date and location of the day's flight 
with a list of significant points keyed to a detailed altitude 
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2 . 


and location plot for the flight. 

A detailed plot of altitudes and locations as a function of 
time after takeoff, with significant points indicated on the plot. 

3. A plot of the extinction, calculated for x = 0.63 pm using the 
assumed optical properties for the particles discussed above for 
the one-minute optical particle data. These and other optical 
data are plotted against time after takeoff for each flight. 

4. A plot of the extinction, calculated for X = 0.63 pm as above 

except for the five-minute optical particle counter data. 

-1 -1 

5. A plot of the backscatter coefficient in m sr 
calculated for A = 0.63 pm based on the measured size 
distribution with five -minute resolution and assumed 
refractive indices. 

6. A plot of the ratio of the calculated backscattering 
coefficient to the calculated extinction coefficient 
for X = 0.63 pm. 

7. A plot of the extinction coefficient calculated for X = 10.6 
pm using the measured size distribution (one-minute averages) 
and assumed refractive indices. 

8. A plot of the extinction coefficient calculated for X = 

10.6 pm as above for the five-minute averages. 

9. A plot of the backscatter coefficient calculated for X = 10.6 
pm using the measured size distribution (five-minute averages) 
and assumed refractive indices. 

10. A plot of the ratio of the calculated backscattering 
coefficient to the calculated extinction coefficient 
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for A = 10.6 jam. 

In addition, a numerical listing of the dN values (in 
cts/cm) for each size interval, presented as 20 minute averages, 
is shown in Appendix B. Although the complete series of plots is 
shown in Appendix A, we will discuss some points of interest in 
the data. 

A representative comparison of the total distribution 
volume, calculated in a manner analogous to that of g and o , 
with the extinction is shown in Figures 3.3 and 3.4. This 
comparison indicates that the variation in particle volume is 
quite comparable to the variation for a £ . Although no formal 
statistical tests have been made, the variation within flight 
legs that is shown in the figure appears to be separable into a 
random variation that is a statistical fluctuation due to the 
small number of particles sampled and a larger scale variation 
that reflects differences in aerosol concentrations on time sales 
of the order of 30 minutes, which correspond to distance scales 
of 200-300 km. There is also, of course, the large variation in 
concentration that is seen in going from the troposphere to the 
boundary layer; this concentration variation results in different 
calculated optical properties for the free troposphere and the 
boundary layer. 

The calculated 3 values also show a variation that is 
similar to that shown by the a £ values; a consideration of 
extinction to backscatter ratios for these distributions, such as 
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EXT. 


***Kno I I enbero Goto*** 
DflTE=MRY 2.1978 



TIME (HOURS AFTER TAKE OFF) 


Figure 3.4. A plot of particulate extinction calculated for 
X = 0.63 pm for the May 2, 1978 data set. 
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the calculations shown in Figure 3.5, shows much less variation 
in calculated ratios. This relative uniformity is due to the 
similar variation of each of these quantities with particle size, 
and to the importance of the submicron particles in all of our 
measured distributions. The major difference that is seen is 
associated with the transition between boundary layer and free 
troposphere. The 10.6 pm values of 3 and o £ are generally 
significantly lower than the 0.63 values, as expected. 
Calculated 3 / a £ ratios for 10.6 pm are near 0.01 and are 
generally somewhat less than those calculated for visible 
wavelength data. 

3.6 LATITUDINAL PROFILES ALONG AIRCRAFT FLIGHT TRACKS 

The X = 0.63 pm extinction data from the individual flights 
were grouped into boundary layer and free tropospheric data and 
plotted against latitude for each of the flight series (1977 and 
1978). The 1978 data were further divided into marine and 
continental data. We note, however, that the separation into 
boundary layer and free troposphere was not always as well 
defined over the continental areas as over the ocean areas. 

The boundary layer data is shown in Figures 3.6 (1977 data), 

3.7 (1978 oceanic data), and 3.8 (1978 continental data). The 
1977 boundary layer data consists entirely of oceanic data. for 
these figures, the symbol 0 represents data between 100 and 200 
m, A data between 200 and 300 m, + data between 300 and 400 m, * 
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RUT ;0<BK5C/SIGMR) 


***Kno I I enberg Data*** 
DRTE=MflY 2.1978 



Figure 3.5. A plot of the ratio of backscatter to extinction 
at 0.63 jjm for the May 2, 1978 data set. 
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LATITUDE (DEG. ) 


Figure 3. 


. Latitudinal profile along the GAMETAG flight 

tracks for the 0.63 pm extinction determined in the 
marine boundary layer during the 1977 GAMETAG 






Extinction 


***Kno i I enberg Qoto*** 
OCEANIC DATA 1978 



Lotitude (Deg.) 


Figure 3.7. As in Fig. 3.6, except for the 1978 oceanic data. 
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data between 400 and 500 m, and Q data between 500 and 1000 m. 
The 1977 and the 1978 data are consistent; each suggests a marine 

-5 -1 

boundary layer extinction somewhat greater than 10 m at 0.63 
pm, with slight decrease south of ITCZ (which was about 5-10 N 
for both of the flight series. There are no clear differences in 
the northern and southern hemisphere data, there are, however, a 
few data points in which the extinction is very low, ~ 2-3 x 
-6 -1 

10 m , which corresponds to very low concentrations of sea salt 
aerosol . 

The free tropospheric data is shown in Figures 3.9, 3.10, 

and 3.11. The 1977 free tropospheric data is shown in Fig. 3.9. 

o 

For this data set, the data for latitudes north of ~ 35 

represent data collected over continental areas; the data for 

o 

latitudes south of 35 are oceanic data. The 1978 free 

tropospheric data is separated into continental data (Fig. 3.10) 
and marine data (Fig. 3.11). The lack of data for latitudes north 
o 

of 20 N in Fig. 3.11 corresponds to missing data between San 

o 

Francisco and Hawaii (20-30 N) which a separate analysis has 
shown to be similar to that measured in 1977. 

For this free tropospheric data, the A's represent data at 
all altitudes between 2 and 3 km, the + 's between 3 and 4 km, the 
*'s between 4 and 5 km, and the Q's for altitudes greater than 5 
km. 
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LATITUDE (DEG.) 


Figure 3.9. Latitudinal profile along the GAMETAG flight 

tracks in the free troposphere for the X =0.63 pm 
extinction data determined during the 1977 flights. 
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9, except for 1978 continental data 





00 15 30 45 60 75 
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for the 1978 oceanic data 





A comparison of the free tropospheric data for the 1977 and 
the 1978 flights indicates that the continental data for each of 
the two years are consistent, with each exhibiting a similar 

-7 

(though wide) range of observed values between ~2 x 10 and 1 x 
-4 -1 

10 m . Both of the marine free tropospheric data sets show a 
decrease in extinction from north to south. Looked at in detail, 
however, there are some significant differences in the two data 
sets. in the latitude region where there is considerable overlap 
o 

of data (25 S to the equator), there appears to be a greater 
extinction in the fall of 1977 than in the spring of 1978 by a 
factor of ~ 2. The 1977 data also show more very high values of 
extinction. While some of these north of the equator may be 
influenced by clouds, others are associated only with increases 
in the submicron aerosol population. The marine data north of 
the equator appear to be consistent for the two years. 

The lowest values of extinction measured in the southern 
hemisphere are lower than the corresponding northern hemisphere 

o 

low values. The 1978 data for latitudes south of 20 S suggest 
that we are measuring a mid-tropospheric background aerosol . The 

o -7 -1 

average for a e for latitudes south of 20 S is ~ 1 x 10 m .If 
we look at the variation about this average, we see that the 
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lowest values measured in the southern hemisphere are 6-7 x 
-8 -1 

10 m with the 200 - 300 km spatial integration of Fig. 3.11, 

-8 -1 -9 -1 

decreasing to 4 x 10 m for the 50 km data set and 4 x 10 m 

for the 9 km data sets. 

As an additional comparison, we have calculated similar 
latitudinal averages of a £ for the 1977 data set for A = 10.6 

pm. These values, shown in Fig. 3.12, are significantly lower 
than the visible wavelength extinction data. The continental 

-8 -1 -5 -1 

10.6 pm data range between 1 x 10 m and 10 m with the 

-6 -1 

majority less than 10 m . The oceanic data shown on this figure 

-9 -1 -5 -1 

ranges from 4x10 m to 10 m , although almost all of the 

-8 -1 

oceanic values are less than 6 x 10 m The average for the 

o 

oceanic data south of about 15 N appears to be in the range of 
-8 -1 

10 m . Because of the greater dependence of the 10.6 pm data on 
the larger particles, and the greater variability of these 
particles, there is more variation in this IR data than in the 
visible data. Based on the data, however, we can estimate that 
the 10.6 pm extinction is lower than the 0.63 pm extinction by a 
factor of about 30. For the more southerly latitudes measured in 

-9 -1 

1978, the 10.6 pm mid-tropospheric extinction is -10 m with 
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***KNOLLENBERG DATA 1977*** 



LRTITUDE (DEG.) 


Figure 3.12. Free tropospheric latitudinal profiles along GAMETAC 
flight tracks for the X = 10.6 pm extinction data. 
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-10 -1 

lowest values (50 km averages) of the order of 3 x 10 m . The 

-7 -1 

corresponding 0.63 urn extinction is ~ 1 x 10 m , roughly 2 
orders of magnitude higher. For the average free tropospheric 
$/CT £ ratio of less than 0.01 calculated for the 10.6 pm data, a 

-11 -1 

average 3 's will be approximately 10 m and lower. 

3.7 CONCLUSIONS 

We have used measured particle size information and assumed 
refractive indices to calculate the aerosol extinction and 
backscatter coefficients along aircraft flight tracks for A = 
0.63 pm and A = 10.6 pm. These calculations show that 

mid-tropospheric extinctions over the Pacific are significantly 
higher in the northern hemisphere than the southern hemisphere. 
The data indicate that the average A = 0.63 pm extinction is - 1 
-7 -1 

x 10 m in the southern hemisphere at the mid-tropospheric 

flight altitudes. The 10.6 pm extinction is between 1 and 2 

r 

tip 

-9 

orders of magnitude less, decreasing to 1 x 10 at the 

southernmost latitudes of the GAMETAG measurements. 
Corresponding calculations of g show that the ratio of 3 to a £ is 
in the range of 0.01 to 0.1 for the 0.63 pm data and in the range 
of 0.01 and lower for the 10.6 pm data. The ratios appear to 
decrease at southern latitudes over the Pacific due to the 
decreased importance of particles with r > 1 pm. 


j 


3-29 



4. MICROPHYSICAL PROCESSES IN AN AEROSOL PLUME 


4.1 BACKGROUND 

A moving aerosol plume may spontaneously undergo different 
microphysical processes. These microphysical processes, in 
addition to condensation and evaporation, include coagulation, 
sedimentation, and diffusion. In our previous work we examined 
the possible effects of growth and evaporative processes on the 
backscatter of radiation traversing an aerosol medium (Deepak et 
al., 1982). In this report, the analysis is focused on the 
effects of coagulation, sedimentation, and diffusion on the 
aerosol size distribution and backscattering coefficient at 10.6 
pm wavelength. In general, a discussion of these microphysical 
processes in an aerosol plume requires a time dependent 3-D 
model. Detailed modeling like this is beyond the scope of this 
analysis. Instead, we will use a 1-D model in this study. This 
simplification allows us to examine these microphysical processes 
in considerable detail. 

4.2 NUMERICAL MODEL 

To achieve the objective of this investigation, a ten-layer 
model has been developed for study of the effect of coagulation, 
sedimentation, and diffusion processes on aerosol size 
distribution and backscattering at 10.6 pm wavelength. Figure 
4.1 shows schematically the model layers. It includes ten 1 km 
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Figure 4.1. Schematic diagram shows the 10-layer model used 
for numerical study of the aerosol transport. The 
arrows indicate the interaction of aerosol size dis- 
tribution between successive layers through the 
effect of gravitational sedimentation. 
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thick layers above the earth's surface. A description of the 
basic equations governing the microphysical processes are given 
below. 


4.2.1 Coagulation 

Coagulation is the process under which aerosol particles 
come into contaqt and coalesce or adhere to one another. In 
other words, coagulation leads to the increase in the number of 
large particles at the expense of the number of smaller 
particles. Coagulation can arise as a result of Brownian motion 
of the aerosol particles — the so-called thermal coagulation. In 
addition, it can be also produced by hydrodynamic, electrical, 
gravitational, or other forces (Fuchs, 1964)* In this 
investigation, we shall assume that the Brownian motion is the 
major cause of the coagulation processes and only the thermal 
coagulation is considered in the analysis. A brief theory of 
this coagulation is given as follows: 


The coagulation equation can be written as 

»oo 


_ 9 _ 

at 


n (V, t ) = 



K(V,u) n(V,t) n(u,t)du 
fV 

K(u i ,Uj) n(u i ,t) n (u j , t ) dV 


0 


(4.1) 


where K(V,u) is the coagulation kernel, which is the frequency of 
collisions per unit volume between particles of volume V and 
particles of volume u, n (V,t) is defined so that the number of 
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particles per unit volume with volume between u and u+du at an 
instant t is given by 
u+du 

J n(u, t ) du. 
u 

The first term on the right-hand side of Eq. 4.1 describes the 
reduction in number of particles with volume V by coagulation 
processes between particles of this particular size with all 
particles of other sizes. The second term on the right-hand side 
of Eq. 4.1 gives the production of new particles with volume V by 
coagulation processes between particles with volumes u and u , 

i j 

so that u and u are related by the equation 

i j 

u + u = V. 

i j 

To calculate the coagulation kernel K(V,u), we have employed the 
formulation derived by Fuchs (1964), whiah is given by 


K'Ui'Uj) - 4* r ij D ij T-tl- + 


4D ii 

G ij r ij 


(4.2) 


with r , D , 6 , and G defined as 

ij ij i j ij 

r = r + r 

ij i j 


D = D + D 

i j i j 

2 2 1/2 
G = (G + G . ) 

ij i 1 
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2 2 i/2 

6 “ ( 6 + 6 ) 
ij i j 

where r and r are the radii of ith and jth particles with 

i j 

volume u and u respectively; D is the diffusion coefficient, 
i j i 

and can be determined from the Einstein relation. That is, D = 

i 

kTB , where k is the Boltzmann constant, T is the absolute 
i 


temperature, and B is the mobility that can be calculated by 

i 


B. = 


i 6 tt i) r 
m 


1 + 1.246 K + 0.24 K exp (- - v 
n n c K 


&] 

n J 


where n is the viscosity of air; and K is the Knudson number 
m n 

and is defined as the ratio of the effective mean free path of 

air molecules {%) to the radius r . G is the average kinetic 

i i 


velocity of a particle, 

1/2 


G . 

l 


f SkT ' 


where m is the mass of the particle of volume u. 


6 . 

x 


is a 


correction factor given by 
1 


^i 6r.il 


• A; i 

i b 


(2r ± + il b ) 3 - (4r ± 2 + Jl 2 ) 


3/2 


2r 


where A = 8D ± / ( 7 rG ± ) 


The analytic solution of Eq. 4.1 is very difficult. In this 
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modelling analysis, it is solved numerically, 
given as follows. 


The approach is 


First, the integration on the right-hand side of Eq. 4.1 is 
achieved by binning the particle volume spectrum into discrete 
bins. In this analysis, the bins are chosen so that the volume 
of a bin is twice that of the preceding bin. To assure the 
conservation of mass as coagulation takes place, it is assumed 
that , 

u i + Uj = f i3 u ± + (1 - f i;1 )u i+1 , jM (4.3) 

Note, f = o, in the case i = j. 
i j 

Equation 4.1 can now be rewritten in the form 


AN i 

At 


= - l K i 3 N i N j 


2 K ij N i 


2 


+ T f ■ K . -N -N . 
j<i ^ 13 i 1 


+ i k . , , . 

2 l-l, l-l 


N 




i N . i N. 


JK, 


(4.4) 


By applying Eq. 4.4, the rate of change of aerosol concentration 
with time in each bin can be calculated, and the change in the 
aerosol size distribution with time can than be determined. In 
the model , the volume range corresponds to a radius range 
of .01 pm to 1.29 pm. The number of bins used in the calculation 
is 32 . 
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4.2.2 Theory of Differential Settling - Sedimentation 
Assume the aerosol particles are spherical in shape. The 


terminal velocity v^. is governed by the Stokes law 


\ = ! t (P E 


r 

ui 


(4.5) 


where 


and 


are densities of the particle and air. 


p m 

respectively, and r is the radius of the particle. Since the 
terminal velocity is a function of particle size, sedimentation 
results in differential settling, and thus in changes of the 
aerosol particle size distribution. Furthermore, the settling of 
the particles leads to interaction between different levels. In 
the present analysis, the fraction of the total number of the 
same size of particles falling out of a given layer is determined 
by 

f = At (4.6) 

r t 

where x is the mean time required for a particle of a given size 

to fall out of the layer and is given by 

t _ _1_ (4.7) 

V t 

and At is the time interval used in the numerical integration. 


3 

In the model computation, we have specified P = 1.001 gm/cm , 


and used U.S. Standard Atmosphere (1976) as the background 
atmosphere. Table 4.1 gives the terminal velocity as a function 
of the particle size. 


4.2.3 Vertical Diffusion 


In one-dimensional models, the vertical movement of 
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Table 4.1. Terminal Velocity V t 



y(|jm) 

V t (in/sec) 

1 

1.00E-02 

1.48E-06 

2 

1.26E-02 

2.35E-06 

3 

1.59E-02 

3.72-E-06 

4 

2.00E-02 

5.91E-06 

5 

2.52E-02 

9.38E-06 

6 

3.17E-02 

1.49E-05 

7 

4.00E-02 

2.36E-05 

8 

5.04E-02 

3.75E-05 

9 

6.35E-02 

5.96E-05 

10 

8.00E-02 

9.45E-05 

11 

1.01E-01 

1.50E-04 

12 

1.27E-01 

2.38E-04 

13 

1.60E-01 

3.78E-04 

14 

2.02E-01 

6.00E-04 

15 

2.54E-01 

9.53E-04 

16 

3.20E-01 

1.51E-03 

17 

4.03E-01 

2.40E-03 

18 

5.08E-01 

3.81E-03 

19 

6.40E-01 

6.05E-03 

20 

8.06E-01 

9.61F-03 

21 

1.02E+00 

1.52E-02 

22 

1.28E+00 

2.42E-02 

23 

1.61E+Q0 

3.84E-02 

24 

2.03E+00 

6.10E-02 

25 

2.56E+00 

9.08E-02 

26 

3.23E+00 

1.54E-01 

27 

4.06E+00 

2. 44E-Q1 

28 

5.12E+00 

3.87E-01 

29 

6.45E+00 

6.15E-01 

30 

8. 13E+00 

9.76E-01 

31 

1.02E+01 

1.55E+00 

32 

1.29E+00 

2.46E+00 
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stratospheric minor constituents is generally described by eddy 
diffusion processes. For small particles this can be an 
important transport mechanism. However, for relative large 
particles the gravitational sedimentation is perhaps the dominant 
transport mode. Generally, the vertical flux at the kth level of 
the aerosol particles due to eddy diffusion can be written as 

= ~ D n rn "a! ( H 2 ] (4,8) 

m 

where n is the concentration, and and denote the air and 

m p 

aerosol, respectively; D is the eddy diffusion coefficient. In 

the model analysis, the set of diffusion coefficients given by 

Liu et al . (1984) is adopted for the computation. They derived 

the vertical diffusion coefficients for tropospheric application 

based on the results of Radon 222 tracer observations. They also 

showed the seasonal dependence of the diffusion coefficients. 

With the vertical flux given by Eq. (4.8), the increase of the 

aerosol particles with size u in the layer between the kth and 

i 

(k -l)th levels in a time step .a t can then simply be calculated 
by 

An i = At (4.9) 

In the calculation, the profile of n is adopted from the 

m 

U.S. Standard Atmosphere (1976). 

4.2.4 Backscattering Coefficient 

The general expression for the volume backscattering 
coefficient p of aerosol particles illuminated with light of 
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wavelength X is 


g (180°) 


oo 

a (X , r,m, 180° ) f(r)dr 

0 


(4.10) 


where r is the radius of the particle, m = (n-ik) the 
refractive index, f(r)dr the number of particles per unit 

o 


with radius between r and r+dr, and o( X, r, m, 180 ) 

backscattering cross section 

x 2 i 1 (X,r,m,180°) + i 2 (X ,r ,m, 180° ) 



complex 

volume 

is the 
(4.11) 


where i and i are the Mie intensity distribution functions for 
1 2 

o 

light scattered at 180 . In the model, the integration over the 
particle size in Eq. 4.10 is carried out numerically. In doing 
this, the computated size spectrum is first fitted by a single 
log-normal or a bimodel aerosol distribution depending on the 
size spectrum. Then, the backscattering coefficient is 
determined with the size parameters obtained from the fitting 
processes. It will be shown latter that a bimodel aerosol 
distribution is a better description in many cases, especially 
when the mixing of aerosol particles of two quite different size 
distributions is taking place. 


4.3 NUMERICAL RESULTS AND DISCUSSION 


The model described in the preceding section has been used 
to examine the evolution of aerosol size distribution and the 
backscattering of an aerosol plume involving coagulation. 
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sedimentation, and diffusion processes. The plume is assumed to 
extend from an altitude of 2 km to 4 km. In the regions below 1 
km and between 4 km and 10 km, only the background aerosol 
particles are introduced in the model . The initial aerosol size 
distributions in various layers are given in terms of parameters 
of log-normal size distributions (Table 4.2). 

th 

In the analysis, it is assumed that the top layer (10 

layer) is steady. In other words, the aerosol size distribution 

in the top layer does not change with time. In the calculation, 
a time step of 0.05 day is used. The solution from this time 
step has been compared with that of a time step of 0.1 day. It 
is found that the differences between the results obtained with 
At = 0.1 day and that with At = 0.05 day are less than 0.5% in 

all the layers considered. The discussion on the computed 
results of the first ten time steps using a 0.05 day increment 
are given below. 

The time evolution of the aerosol size distribution in the 
first layer is given in Figure 4.2. In the figure, the initial 
size distribution is given by the solid line. Only the results 
of every other time step computations are shown. It is 
interesting to note that the concentration of the aerosol 

particles with size greater than about 0.05 pm show a sequential 

increase during the calculated ten steps. On the other hand, a 
decrease in the concentration of particles with size less than 
0.05 pm is also noticeable. This decrease in the concentration 
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Table 4.2. Initial Aerosol Size Distribution 


No. 


(particles/cm ) , 


r 


g 


(mode radius). 


a (width of log-normal curve) 


Layer No. 

Altitude (km) 

No. 

r g l^m) 

0 

1 

0 - 1 

7.0 

0.3 

2.512 

2 

1 - 2 

10.0 

1.0 

2.0 

3 

2-3 

10.0 

1.0 

2.0 

4 

3 - 4 

10.0 

1.0 

2.0 

5 

4-5 

3.1 

0.151 

1.5 

6 

5 - 6 

3. 1 

0.151 

1.5 

7 

6~7 

3.1 

0.151 

1.5 

8 

CO 

t 

3.1 

0.151 

1.5 

9 

8 - 9 

3.1 

0.151 

1.5 

10 

9-10 

3.1 

0.151 

1.5 
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of smaller size particles can be attributed to coagulation, 
sedimentation, and diffusion processes. As to the increase in 
the number concentration of larger particles in the first layer, 
sedimentation and diffusion must be the main mechanisms since 
there is a dust layer right above this first model layer. The 
time evolution of the size distribution of the second layer is 
shown in Fig. 4.3. The main feature of this figure is the 
increase in number concentration of particles with radius smaller 
than 0.05 pm. Very little change is found for larger particles. 
The increase in concentration of smaller particles is mainly due 
to diffusion processes which bring in these smaller particles 
from the layer immediately below. The corresponding time 
evolution of the size distributions of layer No. 3 is given in 
Fig. 4.4. It shows that little change has occurred in the first 
ten time steps. However, a slightly increase in the number 
concentration of small particles is noticeable. Fig. 4.5 shows 
the time variations of the size distributions for the layer No. 
4, which is the upper part of the entire dust layer of this model 


computation . 

The 

noticeable 

features 

of 

Fig. 

4 . 5 are the 

enrichment of 

smal 1 

particles 

(< 0.03 

pm) , 

and 

reduction of 


larger particles. These features can only arise through 
sedimentation and diffusion processes; not aerosol coagulation 
processes, since the latter will cause the changes in opposite 
directions. By inspection, one may find that the increase in the 
number of small particle in the layer No. 4 can be linked with 
the background aerosol particles in the layer No. 5 (Fig. 4.6). 
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Fig. 4.6 displays the time variation of the aerosol size 
distribution in the layer No. 5. The increase in large size 
particles (r > 0.04 pm) are clearly shown. By inspection of the 
aerosol size distributions in the layer No. 4, it is clear that 
the appearance of the larger size particles is primarily due to 
diffusion processes which bring in these particles from the model 
dust layer. From Fig. 4.6, a slight reduction in the number of 
background aerosol particles is also noticeable, which can be 
linked to the increase of the number concentration of particles 
of the corresponding size in the layer No. 4. Since, like the 
layer No. 5, the aerosol particles in the layers above No. 5 are 
specified as background aerosols with the same size distribution, 
the behavior of their aerosol size distributions (not shown) are 
found to have similar features to those of layer No. 5, but of 
much less intensity. It should be mentioned that, as clearly 
shown in Figs. 4.S and 4.6, there is mixing of aerosol particles 
with two quite different size distributions taking place in the 
layers No. 4 and 5 as a result of sedimentation and diffusion 
processes. For this reason, a bimodel fitting process is 
introduced in the model computation. This is done by fitting the 
aerosol size spectrum with a single log-normal mode first. Then 
the differences in the differential aerosol concentrations 
between the original size spectrum and the fitted curve at the 
assigned particle radii (Table 4.1) was determined. A scan 
procedure is then applied to search the maximum difference. A 
second curve fitting process is performed if this maximum 
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Table 4.3. Results of Curve Fitting at the Tenth Time Step 
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-3 -3 

difference is greater than 5 x 10 (cm ). Table 4.3 presents the 
results of this fitting processes after ten steps of 

computation. Figs. 4.7 to 4.11 are the plots of the fitted size 

distributions for model layers from No. 1 to No. 5, 

respectively. By comparing them with the corresponding Figs. 4.2 
to 4.6, one can see that the fitted size distributions correspond 
closely to the original size distributions. 

As mentioned earlier, (Sec. 4.1), one of the objectives of 

this study is to examine the possible effects of the coagulation, 

sedimentation, and diffusion processes on the aerosol 

backscattering at CO wavelengths (10.6 pm). The time evolution 

2 

of the vertical profile of the aerosol backscattering 
coefficients is given in Fig. 4.12. In Fig. 4.12, the symbols O, 
□ r ^ ' k , and Di correspond to the computed results at time 
steps of 0, 2, 4, 6, 8, and 10, respectively. It shows high 

values of the coefficients in the layers No. 2 to No. 4, 

corresponding to the model dust layer. In this dust layer, Fig. 
4.12 shows no noticeable changes in the backscattering 

coefficients. On the other hand, changes can be found in the 

layers immediately below and above this dust layer. It is 

interesting to note that the layer No. 6 does not show any change 
until after the 4th time step of the computation. In addition, 
the changes taking place in layers No. 5 and No. 6 are mainly due 
to the vertical diffusion processes which bring up the large 
particles form the model dust layer. 
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Figure 4.8. As Fig. 4.7, but corresponding Figure 4.9. As Fig. 4.7, but corresponding to 
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BACKSCflTTERING COEFF. (KM-1) 


Figure 4.12. Time evolution of the 10.6 jum backscattering 
function in the model analysis. 

O initial profile, 
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4.4 DISCUSSION AND REMARKS 


The aerosol plume model described in this chapter is able to 

simulate the essential effect of the coagulation, sedimentation, 

and vertical transport, on the aerosol size distributions and the 

aerosol backscattering at CO wavelength (10.6 pm). Although the 

2 

model does not include the aerosol growth effect of condensation 
and evaporation processes (Deepak et al., 1982), an attempt has 
been made to apply the developed model to an observational case 
study. A transport of Asian desert aerosol during the spring of 
1979 has been reported by Shaw (1980), based on observations over 
the Hawaiian Islands and trajectory analysis. A dust layer of 1 
km thick was observed with the layer centered at 3.5 - 4 km 

altitude. SAGE sampling (sunrise measurements) swept gradually 

o o 

over the Pacific Ocean from latitude about 12 to 44 N during the 
period from April 21 to 28, 1979 which coincides with the 
observed Asian dust event. SAGE profiles during this period have 
then been examined at altitudes 4 and 6 km. Unfortunately, it 
was found that many of the SAGE profiles stop at higher 
altitudes. As a result, there is not enough aerosol information 
to give an appropriate picture of this particular dust event. 
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5. CONVERSION OF SAGE i/SAM II EXTINCTION TO 10.6 pM BACKSCATTER 


5.1 INTRODUCTION 

Conversion of a 1 pm SAGE I/SAM II extinction cross section 
to the equivalent backscatter function at 10.6 pm requires 
knowledge or assumptions concerning the aerosol size distribution 
and composition, or alternatively of their general optical 
behavior. In an earlier report, Deepak et al (1982) modeled the 
ratio 6 /a^ Q g as a function of aerosol composition for 
log-normal model of varying mode radii. A summary diagram, 
showing the results of this modeling is reproduced in Fig. 5.1. 
The range for the conversion factor 6 /a-^ 00 shown here is 

three orders of magnitude, sufficiently wide to make it of little 
practical value. This range may be reduced somewhat if we 
restrict ourselves to the free troposphere and take account of 

the fact that the bulk of the smaller aerosols are water soluble 

in composition and that the larger aerosols consist mainly of 
dust. We may then suppose 8 /a. to lie between 10 -4 and 

iU • D X • UU 

—2 ^-1 

10 sr . To improve on this range of values, which is still 

too great, we have chosen to develop a bimodal model to describe 

the aerosol . This has been combined with the size distributions 
measured during the GAMETAG experiments and described in Chapter 
3 of this report. 




MODE RADIUS, 


Figure 5.1. Backscatter to extinction ratios g /a 1 qq 

shown as a function of log-normal mode radius. 
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5.2 BIMODAL SIZE DISTRIBUTION AND FITTING TO THE GAMETAG DATA 


Aerosol size distributions measured during the GAMETAG 
flight series have been fitted to a bimodal model of the form 


dN (r) /d log (r) 


A 
+ A 


. exp (-k.Jin (r/r ) 
l I g-^ 

2 

2 exp (-k^n (r/ r g ^) 


(5.1) 


where k^ 


l/(2£n' 2 a_ ) and k~ = 1/ (2 £n 2 a rf ) 
91 * s 2 


In this equation, r is the particle radius and N(r) the 

comulative particle number. a„ , A »r / r ro and o Q are 

1 2 gi g 2 g x y 2 

constants describing the log normal distributions. Computer 
software has been written which will fit Eq. (5.1) exactly to six 
data points which may be in the form dN(r)/d log(r) as a function 
of r or N(r) as a function of r. The program output contains the 
six constants that describe the bimodal distribution. An example 
of its application to simulated data is shown in Fig. 5.2 where 
the experimental points are shown by round dots, the individual 
fitted curves by dashed lines and the total fitted bimodal 
distribution by the solid line. 


The experimental size distributions used are those listed in 
Appendix B and described in Chapter 3 . They represent particle 
counts obtained for twenty-minute flight segments, the nine size 
intervals being as listed in Table B.l, covering the size range 
0.128 - 2.34 pm. Forty-two flight segments were chosen. 
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dN/d(log r) 


10 



Fig.Si2. Example of the application of the bi-modal log-normal 
fitting program to a simulated experimental size 
distribution. 
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measurements being made at altitudes between approximately 5 and 
7 kilometers and representative of the whole range of latitudes 
covered by the two flight series. Data segments where the 
particle counts were obviously enhanced due to the presence of 
haze or cloud were avoided. Three additional size distributions 
were added from earlier published data (Patterson et al., 1980). 
Size distribution data are available for these flight segments at 
a maximum of nine radii, although inspection of the data in 
Appendix B shows that the channels representing the largest 
particle sizes were frequently empty, particularly over the 
southern hemisphere. The remaining seven or eight data points 
are nevertheless still in excess of the number required to solve 
Eq. (5.1). In order to obtain the best estimate of the log-normal 
constants, the size distributions were plotted and the six data 
points (from the possible maximum of nine) most representative of 
the distribution were used as input to the fitting program. In 
almost all cases, an exact fit was either obtained directly or 
after making a small adjustment to the counts in one channel 


only. Two cases were rejected where reasonable fits could not be 
obtained. Figure 5.3 shows histograms of the resultant values of 


a I a 
91 92 



mode, values 


and r 

92 


of 



It can be seen that, for the accumulation 

lie between 1.1 and 2.0 while r lies 

y i 


between 0.06 and 0.20 pm. For the coarse particle mode, the 


distribution is somewhat narrower, most values of o lying 

y 2 

between 1.0 and 1.5 with r values peaking at around 0.6 pm. 

92 

The somewhat greater imprecision for the accumulation mode is at 
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a r (pm) 

^2 9 2 


Figure 5.3. Histograms of the log-normal parameters obtained from the 
bimodal fitting program. 


(a) 

°9 1 * 

accumulation mode 

(b) 

r 9 1 

accumulation mode 

(c) 
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coarse particle mode 
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coarse particle mode 
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least partially due to the fact that the experimental size ranges 
did not cover particle sizes less than 0.126 pm. 


5.3 BIMODAL OPTICAL MODEL AND COMPARISON OF SAGEI/SAMU AND 

GAMETAG DATA 

The bimodal model described in the previous section has been 
used in conjunction with the AGAUS Mie scattering program to 
calculate the optical properties of the GAMETAG size 
distributions. We have assumed the log normal mode of smallest 
radius (accumulation mode) to consist entirely of water soluble 
aerosol and the mode of greatest radius (coarse particle mode) to 
consist entirely of a dust-type aerosol . The refractive indices 
at 1.00 pm and 10.6 pm are those used previously by Deepak et 
al . (1982) and are listed in Table 5.1. This model is very 

similar to, but not identical to, that used by Georgia Institute 

of Technology to calculate the optical data presented in Chapter 
3 and Appendix A. In the latter case, a two-composition model is 
used, the division between the two components being made at a 

particle radius of 0.5 pm rather than in terms of 

independently-fitted log-normal functions. Comparison of optical 
properties calculated by the two alternative models shows good 
agreement, any differences being much less than observed changes 
due to natural aerosol variation. 

In order to compare the GAMETAG and SAGE I/SAM II data sets, 
the forty- three size distributions described in Section 5.2 have 
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Table 5.1. Refractive Indices 


Aerosol 

Material 



Dust-like 


1.520 


Water 1.520 

Soluble 
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been used in conjunction with the above model to calculate 
corresponding equivalent extinction values at 1 pm. The results 
of these calculations are shown in Fig. 5.4(a) and (b) for the 
1977 and 1978 flights, respectively. As the 1977 flights were 
made in May and June and the 1978 flights in August and September 
(September 1st and 2nd only), the data in these figures may be 
compared with that shown in Fig. 2.6(a) (March-May, 1979) and 
Fig. 2.6(b) (June-August, 1979). The 1 pm extinction values 
calculated for the 1977 flight data (Fig. 5.4(a) agree well with 
the direct SAGE I/SAM II measurements for June-August, 1979 (Fig. 
2.6(b). The absolute magnitudes of the extinction values are very 
close and a similar latitude variation occurs in both data sets. 
The 1978 GAMETAG flight data extends over a wider latitude range 
than the 1977 data and shows greater latitudinal variation. 
Among the calculated values for cr northern hemisphere 
values [Fig. 5.4(b)] agree well with the direct SAGE I/SAM II 
measurements [Fig. 2.6(a)]; the southern hemisphere values are 
somewhat lower than the direct measurements by about a factor of 
two. This difference between the two GAMETAG data sets for 
latitudes south of the equator has been noted in Chapter 3. It 
is believed to be due to a genuine variation in the aerosol 
characteristics observed in the two years and does not reflect 
any instrumental uncertainty. 

Bearing in mind the fact that the GAMETAG data was obtained 
during two cross-sectional flights limited in time and longitude. 
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EXTINCTION v ' extinction 



S 


Figure 5.4. 1 fjm extinction values calculated from size distributions measured during 

the GAMETAG flights. 


the agreement between the calculated 1 jam extinction values and 
the SAGE I/SAM II direct measurements must be considered to be 
within the range of expected variations. Based on this 

agreement, we have used conversion factors for & /o ^ Q0 
calculated from the GAMETAG data to correct the 1 jam SAGE I/SAM 
II data to the equivalent backscatter at 10.6 pm. 

5.4 DERIVATION OF A CONVERSION FACTOR AND ITS APPLICATION 

The forty-three bimodal size distributions taken from the 

GAMETAG data set, as described in the previous section, have been 

used to calculate values for the conversion factor 0Q 

The calculated values are shown in the form of a scatter plot as 

a function of o 1 Q0 in Fig. 5.5. The majority of the values lie 

between 10 -4 and 10 -3 sr -1 and there is a small but definite 

trend in the data, larger values of e^°l 00 being associated 

with larger values of 0Q This is to be expected as 

distributions containing relatively more large particles will 

show increased values of both 6^ a l 00 an< ^ °1 00 A regression 

line for g /a (as the dependent variable) on o (as the 

10.6 1.00 1 ' UU 

independent variable), using logrithmic coordinates, has been 

calculated, and is shown in the figure. Seventy-five percent of 

the calculated values lie within a factor of two of this line. 

The data shown in Fig. 5.5 has been obtained at altitudes of 
5-7 km and is strictly applicable only to this altitude range. 

In spite of this, we have chosen to apply it to all SAGE I/SAM II 
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Figure 5.5. Scatter plot of values of 3 10 .6 /a 1.00 as a 

function of ^ 0Q . Data is for an altitude 

of 5-7 km and taken from the 1977-78 GAMETAG 
flights. The figure also shows the best fit 
regression line for g la ^ Q0 upon o ^ 0Q . 
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tropospheric data above about 4 km. Some error may be expected 

at 4 km (and below) because of relatively increased numbers of 

large particles; these will increase 3 /a faster than they 

10.6 1.00 

will increase o At the higher altitudes near the tropopause, 

■i • U U 

the aerosol composition will most likely contain a high 

percentage of sulfuric acid and its optical properties will merge 
into those of the stratospheric aerosol. The size distribution 
will also be strongly dependent upon volcanic input and the 
conversion factor likewise highly variable (Kent et al., 1984). 
Despite these limitations, it is felt that the regression line 
shown in Fig. 5.5 is a good approximation for use, during a 
volcanically quiet period, in the middle and upper free 

troposphere. 

Table 5.2 contains a list of conversion factors 

corresponding to the regression line in Fig. 5.5. This has been 

applied to the contour plots of 1 pm extinction shown in Fig. 

2.8. The results of this are shown in Fig. 5.6(a) and (b) for the 

March-May and September-November , 1979 seasons, respectively. We 

see that, in the middle and upper free troposphere, southern 

-11 

hemisphere values for g range between 3 x 10 and 2 x 

10 In the northern hemisphere, values are greater at all 

altitudes becoming as large as 7 x 10 10 in March-May, 1979 at 5 
km and 70N. The latitudinal and seasonal variation is very 
pronounced, as is the variation with altitude. It must be 
emphasized that these are seasonal median values and daily values 
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Table 5.2. Conversion Factors Used with the SAGE 1/ 
SAM 11 Free Tropospheric Data Sets 


° 1.00 
{id 7 rtf 1 ) 
1.0 


1.5 


2.0 


3.0 


4.0 


5.0 


6.0 


8.0 


10.0 

15.0 


20.0 


3 /a 

10 . 6 / 1.00 

-4 - 1 

(10 sr '} 

1.34 
1.53 
1.67 
1.89 
2.07 
2.22 

2.35 
2.57 
2.75 
3.12 
3.42 
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Figure 5.6. Contour plots of the free tropospheric backscatter function at 10.6 /urn derived from the 
SAGE I /SAM II data set. 


would be expected to depart considerably from them. In addition, 


they correspond to a volcanically quiet period 
present, more volcanically active period, upper 
values would be considerably increased. 


during the 
tropospheric 



6. C0 2 LIDAR MEASUREMENTS OF 
6.1 PUBLISHED MEASUREMENTS 

Published measurements of aerosol backscatter cross sections 

in the free troposphere using CO lidar systems are scarce. The 

2 

most extensive data series is that made at NOAA in Boulder, 
Colorado (Post et al., 1982; Post, 1984a) . Measurements have also 
been made at NASA-MSFC , Huntsville, Alabama (Weaver, 1983; Jones, 
1983) and JPL, Pasadena, California (Menzies et al., 1984). 
Updates to these data sets have recently been presented at the 
Third NASA/NOAA Infrared Lidar Backscatter Workshop, Lake Tahoe, 
Nevada, January 14, 1985. These data sets have all been obtained 
over, or close to, the continental U.S.A. There is, in addition, 
an extensive unpublished data set obtained by the Royal Signals 
and Radar Establishment, United Kingdom over the United Kingdom 
and Europe (Vaughan, 198 5) . , 

The earliest published data that extends to altitudes above 
4-5 km (Post et al., 1982) was taken in the Spring and Summer of 
1981 when the atmosphere was affected by aerosol produced in 
recent volcanic eruptions. The only earlier data, which shows 
conditions prior to these eruptions is that of .Schweisow et al. 
(1981). This data, taken in 1978, does not extend above 5 km (4 
km for most of the data) and thus does not describe conditions in 
the middle and upper free troposphere. Data, since 1981, has 
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been very strongly affected by the eruption of El Chichon in 
April, 1982 and we therefore have no direct measurements of 
background non-volcanic aerosol scattering cross sections in the 
upper free troposphere. 

6.2 VARIATION OF 3, n , WITH ALTITUDE AND SEASON 

Figure 6.1 summarizes the measured variation of g with 

altitude. Figure 6.1(a) shows eight seasonal averages for 8 -^q 6 
measured by NOAA (Post, 1983; Post, 1984a), Figure 6.1(b) shows a 
single seasonal average obtained by JPL (Menzies, 1984) and 
Figure 6.1(c) and (d) show individual profiles measured by the 
NASA-MFC airborne system over California in 1981 and 1982, 
respectively (Jones, 1983). Table 6.1 summarizes the main 
characteristic of the RSRE measurements (Vaughan, 1983 ) . Although 
these data have been obtained by several different instruments 
and at different times and locations, several general features 
may be identified. 

1. Maximum values for ^iq. 6 occur i- n the boundary layer 

— Q 

where , is greater than 10 

1 U • o 

2. Above the boundary layer, g falls rapidly to a rather 

1 U . D 

broad minimum between about 5 and 10 km altitude. Values in this 
minimum lie mainly in the range 3 x 10 -11 - 3 x 10 -1 ® m”* sr -1 . 

3. Above the minimum, values rise into the stratosphere 
where peak values of the order of 10 - ^ itf^ sr"-'- may be found 
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Figure 6.1. Measured values for g. Q g in the free troposphere. 





Table 6.1. Royal Signal and Radar Establishment 

$ nr . Measurements (D. M. Vaughan, 1983) 
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(post-volcanic) . 


Apart from these general characteristics of the vertical 
profile, more specific characteristics of 6 , shown by one or 
more data sets are as follows. 

1. There exists a day to day variation (up to 5 km at least) 
of an order of magnitude or more, combined with pronounced 
vertical stratification. 

2. Geophysical differences exist between California and 
Colorado. Over the former, the decrease in ^ occurs at a 
lower altitude than over the latter. This is most likely 
associated with the higher altitude of the Colorado plateau and 
greater vertical convection. 

3. A seasonal variation is observed by NOAA and by RSRE . At 

an altitude of 6 km minimum values of g, are seen by NOAA in 

Pall and Winter ( ~ 10 10 m 1 sr ^ , maximum values occur in 

-9 -1 -1 

Spring and Summer ( ~ 10 m sr ). RSRE observed minimum g 
in Winter. 

4. A decrease occurred in the peak altitude of the 
stratospheric volcanic layer over NOAA (from 23 km in Fall, 1982 
to 17 km in Fall, 1983). 

5. Lowest mean values for gin the upper stratosphere (2 
x 10 11 m ^) were observed by NOAA in Spring and Summer of 1982. 
This was after the decay of material injected by the eruptions of 



St. Helens and Alaid and before the arrival of material injected 
by El Chichon. 

6.3 PROBABILITY DISTRIBUTIONS FOR , 

10.6 

Cumulative probability distributions for g have been 

published by NOAA (Post et al. f 1982; Post, 1984a). These show a 
good fit to log-normality and the best-fit straight lines to the 
data sets at five altitudes and on two occasions are reproduced 
in Fig. 6.2(a) and (b) (Post et al., 1982; Post, 1984a, b). The 
lower spread in 6 values at stratospheric altitudes, as 
opposed to tropospheric altitudes is shown in Fig. 6.2(b) and is 
to be noted. 
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7. DATA INTERCOMPARISON 


7.1 VARIATION WITH LATITUDE AND SEASON 

Figures 7.1(a) and (b) are composite figures on which 
theoretical and experimental values for 3^0 6 at an altitude of 
about 6 km have been plotted as a function of latitude. The 
altitude of 6 km has been chosen for several reasons. It is 
close to the mean flight altitude for much of the GAMETAG 
measurement series and it is an altitude at which good (> 40%) 
SAGE I/SAM II penetration is found, and for which the median 
values may be considered to be without serious systematic bias. 
It is also close to the level at which minimum values for 3 10 g 
are found experimentally and it thus represents a compromise 
altitude that is not too affected by ground topography below and 
volcanic fallout from above. 

Data for Spring and Summer have been plotted in Fig. 7.1(a), 
the much more limited data for Fall and Winter is plotted in Fig. 
7.1(b). In Fig. 7.1(a), the GAMETAG data points represent median 
values for 3 -^q g observed over flight segments for which the 
flight altitude was between 5 and 7 km. Data for 1977 and 1978 
have been distinguished by different symbols. The lines 
representing the SAGE I data are the March-May and 
September-November median values in Fig. 7.1(a) and 7.1(b), 
respectively . 
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Seasonal averages for JPL and NASA-MSC are shown by single 
symbols (no data is available for Fall-Winter). The larger amount 
of data from NOAA is shown by vertical lines in Figs. 7.1(a) and 
7.1(b) which represent the range of seasonal mean values found 
[five Spring-Summer seasons in Fig. 7.1(a) and three Fall-Winter 
seasons in Fig. 7.1(b)]. 

Examination of these figures shows that, apart from the 1978 

GAMETAG data in the southern hemisphere, agreement is quite good 

between the different data sets. Comment has already been made 

concerning the difference between the 1977 and 1978 GAMETAG data 

sets and without further evidence we must regard this as a 

sampling fluctuation. In the northern hemisphere, most of the 

-10 

measured and modeled values for « lie between 10 and 2 x 

p 10 . 6 

-9 -1 -1 

10 m sr When judging the goodness of fit between these 

data sets, it must be remembered that they refer to different 
years and geographical regions and to different degrees of 
stratospheric volcanic contamination. There is apparently still 
a very considerable need for direct measurements of $20.6 i n t ^ ie 
southern hemisphere. 

7.2 VARIATION WITH ALTITUDE 

Figure 7.2 shows the altitude variation of 6 for three 
early NOAA seasonal averages, the 1983-1984 JPL Spring-Summer 
average and the SAGE I modeled values for 40 N. One of the NOAA 
averages (Summer 1981) and the JPL average show a clear 
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Figure 7.2. Composite diagram showing modeled and 

directly measured vertical profiles for 3 1Q 
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stratospheric volcanic enhancement which is not present in the 
other profiles. Between altitudes of 5 and 9 km, the agreement 
between SAGE I and JPL is good. The agreement between both of 
these and the NOAA data is quantitatively poor. it is quite 
possible that much of this difference may be due to the location 
of Boulder high on the continental plateau. There will be 
increased convective activity and the boundary layer will extend 
to a considerably greater altitude above sea-level. The low 
values observed by NOAA at a height of 9 km in Spring and Summer 
of 1982 do not fit well with the modeled data. 

7.3 PROBABILITY DISTRIBUTIONS 

Examples of the probability distribution for a and p 

-L »UU 1 v • O 

have been shown in Figs. 2.3 (a) -(e) and 6.2(a) and (b). Certain 

similarities exist in these figures. All the NOAA data shows 
clear log-normal characteristics. This is also shown by the SAGE 
I/SAM II data in the stratosphere. In the upper troposphere, the 
distribution is clearly two-component and we believe the higher 
estimation component to be sub-visible cirrus (cirrus returns 
have been omitted from the NOAA data). In the middle free 
troposphere, the SAGE I/SAM II data is once again approximately 
log-normal for some of the latitude bands. In order to make a 
more quantitative comparison of the distribution, Table 7.1 shows 
the extinction and backscatter changes corresponding to the 5% 

50% probability range, i.e., the departure of the 5% probability 
level from the median level. Values are shown for similar 



latitudes (20 N - 60 N for SAGE I and 40 N for NOAA) and for 
altitudes above 5 km. The range of the optical parameter has 
been expressed as the ratio of its value at the 50% probability 
level to its value at the 5% level. 

Table 7.1 shows that for both C]_ qq and g / the ratio 
decreases with altitude by roughly comparable amounts. The 
absolute magnitude of the ratio in the two cases is nevertheless 
very different. Part of this difference can be accounted for by 
the fact that the conversion factor g /o^ 0Q itself depends 
upon the value of © This is, however, not sufficient to 
account for most of the observed difference. It is more likely 
that the explanation for this lies in the different geometry of 
the two measurements. The 10.6 pm backscatter function is 
measured over horizontal dimensions of a few meters only and the 
probability distribution will reflect all changes with horizontal 
scales greater than this. The 1 pm extinction function, on the 
other hand, is derived from a measurement made over a 200-300 km 
long horizontal path. Fluctuations with horizontal scales less 
than this will be smoothed. The major part of the difference in 
ratio observed in Table 7.1 almost certainly arise from these 
different measurement geometries. 

7.4 VOLCANIC EFFECTS 

Volcanic effects have already been discussed in Sections 2.5 
and 6.2 and the details will not be presented again here. The 
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most obvious point of comparison to note is the stratospheric 
enhancement observed in both cases (Figs. 2.n and 6.1). A more 
subtle point of comparison is the enhancement in the upper 
troposphere. Figure 2.11 shows that the upper troposphere, 
following volcanic injection into the stratosphere, has a raised 

I pm extinction which extends from the tropopause down to 5-6 
km. Figure 6.1(a) shows similar effects, post-volcanic 10.6 pm 
back-scattering, following volcanic injection, is enhanced from 
the tropopause down to an altitude of 7-8 km. The fact that this 
altitude is somewhat higher than that given above for SAGE I/SAM 

II is, as discussed earlier, most likely due to a geographic 
effort associated with the location of Boulder, Colorado. 
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8 . CONCLUSIONS 


In the work described in this report, the objective has been 
to improve our understanding of the global behavior of the 
aerosol backscatter function in the free troposphere at a 
wavelength of 10.6 pm. Three experimental data sets have been 
used, each with certain limitations in either wavelength or 
temporal and geographical coverage. In addition to the use of 
these data sets, modeling has been carried out on the 
microphysical and optical changes to be expected within an aerosol 
plume . 

The SAGE I/SAM II data set has yielded a detailed picture of 
the variation of 1 pm aerosol extinction with latitude, season 
and altitude above 5 km. This shows pronounced hemispheric 
asymmetry, with minimum extinction occurring in high southern 
latitudes. In both hemispheres, the aerosol extinction in spring 
and summer is greater than that in fall and winter. Volcanic 
enhancement of the aerosol extinction above 5 km is also apparent 
at high latitudes in the 1980-81 data. The GAMETAG data sets for 
1977 and 1978 have shown a latitude variation in aerosol 
concentration and optical properties, that is in overall 
qualitative agreement with the SAGE I/SAM II data, and in good 
quantitative agreement in the northern hemisphere. Direct 
measurements of 3^o 5 using C0 2 lidar are confined to a limited 
latitude band in the northern hemisphere and have only been made 
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since 1980, during a period of volcanic activity. 


Comparison has made of the direct measurements with values 

for g calculated from the other data sets. This comparison 

has particularly been carried out for an altitude of 6 km, which 

is close to the level at which minimum values for g,„ are 

p 10 . 6 

currently observed. Median summer values for r . at this 

10.6 

level, in the northern hemisphere, vary from 1 x 10 10 m -1 sr - l 

— 9 _ i __ -j 

at the equator to 1 x 10 m sr at 60 N. Median winter values 

at the same altitude are about 2 x 10" at all northern 
latitudes. In the southern hemisphere, the situation is less 
clear. There are no direct measurements and the 1977 and 1978 
GAMETAG measurements yield significantly different values for 
8 1ri ... SAGE I/SAM II median values at 6 km altitude lie between 3 
x 10 - llxn -1 sr -1 and 2 x 10 _10 m~ 1 sr _1 ; the range for median GAMETAG 
values is between 1 x 10 11 and 2 x 10 10 , Both data sets show a 
decreasing aerosol concentration toward high southern latitudes. 
Intercomparison of vertical profiles for g^g g is difficult. 
GAMETAG data has been obtained mainly between 5 and 7 km altitude 
and it is not possible to derive a vertical profile with any 
precision. Almost all the direct measurements of g^ Q g are 
affected by volcanic input, which enhances the upper free 
troposphere. Reasonable agreement is obtained between the direct 
measurements of 6 an< ^ va l ues modeled from SAGE I/SAM II at 

altitudes between 5 and 10 km but any exact comparison is not at 
present possible. 
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It is clear, from the analysis carried out and the data 
comparisons made, that we now have a reasonable quantitative 
understanding of the principal features of free tropospheric 
aerosol concentrations and optical characteristics and their 
variation with latitude season and volcanic activity. We are, 
nevertheless, still a considerable distance from being able to 
supply a complete statistical description, particularly of the 
optical characteristics at the longer wavelengths. This is 
especially true in the southern hemisphere where there is a total 
lack of direct measurements at 10.6 pm. 
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Fig. A 1 (a). GAMETAG flight data for August 7, 1977. 

Altitude and location flight track plotted as a 
function of time after takeoff. 
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Fig. A 1 (f). GAMETAG flight data for August 7, 1977. 

Calculated particulate extinction along the flight 
track for one-minute data sets for X= 10.6 urn. 
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Fig. Al (i). GAMETAG flight data for August 7, 1977. 

Calculated ratios for backscatter to extinction for 
five-minute data sets for X= 10.6 lin. 
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Table A2. Significant times for August 8, 1977. 

Portland, Oregon to Anchorage, Alaska. 
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Fig. A 2 (a). GAMETAG flight data for August 8, 1977. 

Altitude and location flight track plotted as a 
function of time after takeoff. 
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Fig. A 2 (b). GAMETAG flight data for August 8, 1977. 

Calculated particulate extinction along the flight 
track for one-minute data sets for X= 0.63 ym. 
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Fig. A 2 (c). CAMETAG flight data for August 8, 1977. 

Calculated particulate extinction along the flight 
track for five-minute data sets for X= 0.63 ym. 


i 


A14 




BKSC . X l.E-7 (1/M . STR) 


***KNOLLENBERG DRTR*** 
DRTE= 8/ 8/1977 



TIME (HOURS RFTER TRKE OFF) 


Fig. A 2 (d). GAMETAG flight data for August 8, 1977. 

Calculated backscatter coefficient along the flight 
track for five-minute data sets for X = 0.63 yin. 




A15 


mid 




RRTIOc 3K5C/S1GMR) 


***KNOLLENBERG DRTR*** 
DRTE= 8/ 8/1977 



TIME (HOURS AFTER TAKE OFF) 


Fig. A 2 (e). GAMETAG flight data for August 8, 1977. 

Calculated ratios for backscatter to extinction for 
five-minute data sets for A= 0.63 jam. 
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Fig. A 2 (f). GAMETAG flight data for August 8, 1977. 

Calculated particulate extinction along the flight 
track for one-minute data sets for A= 10.6 um. 
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Fig. A 2 (g). GAMETAG flight data for August 8, 1977. 

Calculated particulate extinction along the flight 
track for five-minute data sets for A= 10.6 ym. 
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Fig. A 2 (h). GAMETAG flight data for August 8, 1977. 

Calculated backscatter coefficient along the flight 
track for five-minute data sets for A= 10.6 ym. 
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Fig. A 2 (i). GAMETAG flight data for August 8, 1977. 

Calculated ratios for backscatter to extinction for 
five-minute data sets for A= 10.6 v* 1 ** 
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Table A3. Significant times for August 22, 1977 
Denver, Colorado to Moffett Field, 
California. 

-* %■ 

Significant 

Points 


# 

TIME 


1 

18:22 

Denver 

2 

19:12 


3 

19:22 


4 

20:21 


5 

21:01 


6 

21:15 


7 

21:17 


8 

21:38 


9 

21:41 


10 

21:51 


11 

22:17 


12 

22:26 


13 

23:15 


14 

23:34 

Moffett Field 
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Fig. A3 (a). GAMETAG flight data for August 22, 1977. 

Altitude and location flight track plotted as 
function of time after takeoff. 
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Fig. A3 (b). GAMETAC flight data for August 22, 1977. 

Calculated particulate extinction along the flight 
track for one-minute data sets for A= 0.63 urn. 
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Fig. A3 (c). GAMETAG flight data for August 22, 1977. 

Calculated particulate extinction along the flight 
track for five-minute data sets for A= 0.63 ym. 
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Fig. A3 (d). GAMETAG flight data for August 22, 1977. 

Calculated backscatter coefficient along the flight 
track for five-minute data sets for X= 0.63 ym. 
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Fig. A3 (e). GAMETAG flight data for August 22, 1977. 

Calculated ratios for backscatter to extinction for 
five-minute data sets for A= 0.63 ym. 
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Fig. A3 (f). GAMETAG flight data for August 22, 1977. 

Calculated particulate extinction along the flight 
track for one-minute data sets for X= 10.6 pm. 
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Fig. A3 (g). GAMETAG flight data for August 22, 1977. 

Calculated particulate extinction along the flight 
track for five-minute data sets for A= 10.6 ym. 
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Fig. A3 (h). GAMETAG flight data for August 22, 1977. 

Calculated backscatter coefficient along the flight 
track for five-minute data sets for X= 10.6 ym. 
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Fig. A3 (i). GAMETAG flight data for August 22, 1977. 

Calculated ratios for backscatter to extinction for 
five-minute data sets for X = 10.6 nn. 


J 


A30 




Table A4. 


Significant times for August 23, 1977. 
Moffett Field, California to Hilo, Hawaii. 


Significant Points 


# TIME 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 


17:45 Moffett Field 
18:20 
20:31 
20:40 
22:13 
22:19 
23:18 
23:28 
00:09 
00:26 
00:46 

00:59 Hilo 
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Fig. A 4 (a). GAMETAG flight data for August 23, 1977. 

Altitude and location flight track plotted as a 
function of time after takeoff. 
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Fig. A 4 (b). GAMETAG flight data for August 23, 1977. 

Calculated particulate extinction along the flight 
track for one-minute data sets for X — 0.63 ym. 
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Fig. A4 (c). GAMETAG flight data for August 23, 1977. 

Calculated particulate extinction along the flight 
track for five-minute data sets for A= 0.63 ym. 
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Fig. A4 (e). GAMETAG flight data for August 23, 1977. 

Calculated ratios for backscatter to extinction for 
five-minute data sets for X = 0.63 ym. 
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***KNOLLENBERG DRTR*** 
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Fig. A 4 (g). GAMETAG flight data for August 23, 1977. 

Calculated particulate extinction along the flight 
track for five-minute data sets for X= 10.6 ym. 
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Fig. A 4 (h). GAMETAG flight data for August 23, 1977. 

Calculated backscatter coefficient along the flight 
track for five-minute data sets for A= 10.6 ym. 
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Fig. A 4 (i). GAMETAG flight data for August 23, 1977. 

Calculated ratios for backscatter to extinction for 
five-minute data sets for A= 10.6 ym. 
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Table A5. 


Significant times for August 25, 1977. 
Hilo, Hawaii to Johnston Atoll. 


Significant Points 


# TIME 


1 

20:51 

2 

21:07 

3 

22:03 

4 

22:19 

5 

22:41 

6 

22:54 

7 

23:09 

8 

23:47 

9 

23:59 

10 

00:32 

11 

01:14 

12 

01:33 

13 

02:36 

14 

02:48 


Hilo 


Johnston Atoll 
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Fig. A 5 (a). GAMETAG flight data for August 25, 1977. 

Altitude and location flight track plotted as a 
function of time after takeoff. 
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Fig. A 5 (b). CAMETAG flight data for August 25, 1977. 

Calculated particulate extinction along the flight 
track for one~minute data sets for X = 0.63 pi. 
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Fig. A 5 (c). GAMETAG flight data for August 25, 1977. 

Calculated particulate extinction along the flight 
track for five-minute data sets for X = 0.63 ym. 
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Fig. A 5 (d). GAMETAG flight data for August 25, 1977. 

Calculated backscatter coefficient along the flight 
track for five-minute data sets for A= 0.63 pm. 
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Fig. A 5 (e). GAMETAC flight data for August 25, 1977. 

Calculated ratios for backscatter to extinction for 
five-minute data sets for X = 0.63 ym. 
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Fig. A 5 (f). GAMETAG flight data for August 25, 1977. 

Calculated particulate extinction along the flight 
track for one-minute data sets for X = 10.6 nm. 
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***KNOLLENBERG DRTR*** 
DRTE= 8/25/1977 



Fig. A 5 (g). GAMETAG flight data for August 25, 1977. 

Calculated particulate extinction along the flight 
track for five-minute data sets for X= 10.6 jam. 
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***KNOLLENBERG DRTR*** 
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Fig. A 5 (h). GAMETAG flight data for August 25, 1977. 

Calculated backscatter coefficient along the flight 
track for five-minute data sets for A= 10.6 ym. 
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Fig. A 5 (i). GAMETAG flight data for August 25, 1977. 

Calculated ratios for backscatter to extinction for 
five-minute data sets for A= 10.6 jm. 
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Table A6. Significant times for August 26, 1977. 

Johnston Atoll to American Samoa. 


Significant Points 


# 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 
13 


TIME 

20:00 Johnston Atoll 

20:35 
22:37 
23:01 
23:47 
00:06 
00:41 
00:50 
00:52 
02:20 
02:39 
02:48 

02:55 Pago Pago, American Samoa 
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Fig. A 6 (a). GAMETAG flight data for August 26, 1977. 

Altitude and location flight track plotted as a 
function of time after takeoff. 
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Fig. A 6 (b). GAMETAG flight data for August 26, 1977. 

Calculated particulate extinction along the flight 
track for one-minute data sets for A= 0.63 ym. 
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Fig. A 6 (c). GAMETAG flight data for August 26, 1977. 

Calculated particulate extinction along the flight 
track for five-minute data sets for X= 0.63 ym. 
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Fig. A 6 (d). GAMETAG flight data for August 26, 1977. 

Calculated backscatter coefficient along the flight 
track for five-minute data sets for A= 0.63 ym. 


i 


A55 




RRT I n ' BKSC/S I GMR ) 


***KNOLLENBERG DRTR*** 
DRTE= 8/26/1977 



. 0 1 2 3 4 5 6 7 


TIME (HOURS AFTER TAKE OFF) 


Fig. A 6 (e). GAMETAG flight data for August 26, 1977. 

Calculated ratios for backscatter to extinction for 
five-minute data sets for A= 0.63 ym. 
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***KNOLLENBERG DRTR*** 
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Fig. A 6 (f). GAMETAG flight data for August 26, 1977. 

Calculated particulate extinction along the flight 
track for one-minute data sets for X= 10.6 jim. 
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Fig. A 6 (g). GAMETAG flight data for August 26, 1977. 

Calculated particulate extinction along the flight 
track for five-minute data sets for X= 10.6 ym. 
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Fig. A 6 (h). GAMETAG flight data for August 26, 1977. 

Calculated backscatter coefficient along the flight 
track for five-minute data sets for X= 10.6 ym. 
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Fig. A6 ( i) . GAMETAG flight data for August 26, 1977. 

Calculated ratios for backscatter to extinction for 
five-minute data sets for X= 10.6 ym. 
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Table A7. Significant times for August 28, 1977. 

Pago Pago, American Samoa and Return 
(North) . 


Significant Points 


# 

TIME 

1 

21:02 

2 

21:35 

3 

22:30 

4 

23:00 

5 

23:45 

6 

00:18 

7 

01:24 

8 

01:57 

9 

02:18 

10 

02:20 

11 

02:36 

12 

02:43 

13 

02:55 


Pago Pago 


Pago Pago 
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Fig. A 7 (a). GAMETAG flight data for August 28, 1977. 

Altitude and location flight track plotted as 
function of time after takeoff. 



ExT. x 1 .E-7 ( 1/M) 


ORIGINAL 
Of POOR 


PAGE IS 

Quality 


***KNOLLENBER DOT R*** 
DRTE= 8/28/ 977 



T I ML (HOURS RFT:R TAKE OFF) 


Fig. A 7 (b). GAMETAG flight data for August 28, 1977. 

Calculated particulate extinction along the flight 
track for one-minute data sets for X= 0.63 ym. 
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Fig. A 7 (c). GAMETAG flight data for August 28, 1977. 

Calculated particulate extinction along the flight 
track for five-minute data sets for X- 0.63 ym. 
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Fig. A 7 (d). GAMETAG flight data for August 28, 1977. 

Calculated backscatter coefficient along the flight 
track for five-minute data sets for X= 0.63 pro* 
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Fig. A 7 (e). GAMETAG flight data for August 28, 1977. 

Calculated ratios for backscatter to extinction for 
five-minute data sets for A= 0.63 ym. 
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Fig. A 7 (f). GAMETAG flight data for August 28, 1977. 

Calculated particulate extinction along the flight 
track for one-minute data sets for A= 10.6 nm. 
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Fig. A 7 (g). GAMETAG flight data for August 28, 1977. 

Calculated particulate extinction along the flight 
track for five-minute data sets for A= 10.6 ym. 
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Fig. A 7 (h). GAMETAG flight data for August 28, 1977. 

Calculated backscatter coefficient along the flight 
track for five-minute data sets for X = 10.6 ym. 
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Fig. A 7 (i). GAMETAG flight data for August 28, 1977. 

Calculated ratios for backscatter to extinction for 
five-minute data sets for X= 10.6 vm. 
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Table A8. Significant times for August 31, 1977. 

Pago Pago, American Samoa and Return 
(South) . 

Significant 

Points 

! 

_#_ 

TIME 


1 

20:58 

Pago Pago 

2 

21:30 


3 

22:32 


4 

22:38 


5 

23:44 


6 

23:53 


7 

01:14 


8 

01:35 


9 

01:53 


10 

02:07 


11 

02:25 


12 

02:33 


13 

02:39 


14 

02:45 


15 

03:19 


16 

03:24 


17 

03:39 


18 

03:49 

Pago Pago 
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Fig. A 8 (a). GAMETAG flight data for August 31, 1977. 

Altitude and location flight track plotted as 
function of time after takeoff. 
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Fig. A8 (b). GAMETAG flight data for August 31, 1977. 

Calculated particulate extinction along the flight 
track for one-minute data sets for X= 0.63 ym. 
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Fig. A 8 (c). GAMETAG flight data for August 31, 1977. 

Calculated particulate extinction along the flight 
track for five -minute data sets for X= 0.63 ym. 
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Fig. A8 (d). GAMETAC flight data for August 31, 1977. 

Calculated backscatter coefficient along the flight 
track for five-minute data sets for X = 0.63 ym. 
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Fig. A 8 (e). GAMETAG flight data for August 31, 1977. 

Calculated ratios for backscatter to extinction for 
five-minute data sets for X= 0.63 ym. 
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Fig. A8 (f). GAMETAG flight data for August 31, 1977. 

Calculated particulate extinction along the flight 
track for one-minute data sets for X= 10.6 jam. 
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Fig. A 8 (g). GAMETAG flight data for August 31, 1977. 

Calculated particulate extinction along the flight 
track for five-minute data sets for A= 10.6 jam. 
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Fig. A 8 (h). GAMETAG flight data for August 31, 1977. 

Calculated backscatter coefficient along the flight 
track for five-minute data sets for X = 10.6 ym. 
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Fig. A 8 (i). GAMETAG flight data for August 31, 1977. 

Calculated ratios for backscatter to extinction for 
five-minute data sets for X = 10.6 lin. 
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Table 

A9. Significant times for September 1, 1977. 

Pago Pago, American Samoa to Johnston Atoll. 

< *■ 

Significant 

Points 


# 

TIME 


1 

19:57 

Pago Pago 

2 

20:27 


3 

21:57 


4 

22:15 


5 

22:42 


6 

23:05 


7 

00:49 


8 

01:01 


9 

01:19 


10 

01:31 


11 

01:32 


12 

01:41 


13 

02:01 


14 

02:25 


15 

02:40 


16 

02:46 

Johnston Atoll 


i 

A81 

* 



-0*r* .... 



J 


A82 


Fig. A 9 (a). GAMETAG flight data for September 1, 1977. 

Altitude and location flight track plotted as a 
function of time after takeoff. 
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Fig. A 9 (b). GAMETAG flight data for September 1, 1977. 

Calculated particulate extinction along the flight 
track for one-minute data sets for X= 0.63 ym. 
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Fig. A 9 (c). GAMETAG flight data for September 1, 1977. 

Calculated particulate extinction along the flight 
track for five-minute data sets for X= 0.63 ym. 
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Fig. A 9 (d). GAMETAG flight data for September 1, 1977. 

Calculated backscatter coefficient along the flight 
track for five-minute data sets for X = 0.63 pm . 
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Fig. A 9 (e). GAMETAG flight data for September 1, 1977. 

Calculated ratios for backscatter to extinction for 
five-minute data sets for A= 0.63 ym. 
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Fig. A 9 (i). GAMETAG flight 

Calculated ratios 
five-minute data 






Table A10. Significant times for September 2, 1977 
Johnston Atoil to Hilo, Hawaii. 


Significant Points 


# 

TIME 

1 

16.49 

2 

17:13 

3 

18:02 

4 

18:30 

5 

19:37 

6 

19:48 

7 

20:06 

8 

20:16 

9 

20:36 

10 

21:18 

11 

21:32 

12 

21:47 

13 

22:10 

14 

22:27 

15 

22:40 


Johnston Atoll 


Hilo 
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Fig. A 10 (a). GAMETAG flight data for September 2, 1977. 

Altitude and location flight track plotted as a 
function of time after takeoff. 


TIME (HOURS AFTER TAKE OFF) 


Fig. A 10(b). GAMETAG flight data for September 2, 1977. 

Calculated particulate extinction along the flight 
track for one-minute data sets for X= 0.63 ym. 
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Fig. A 10 (c) . GAMETAG flight data for September 2, 1977. 

Calculated particulate extinction along the flight 
track for five-minute data sets for X= 0.63 pm. 
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Fig. AlO(f). GAMETAC flight data for September 2, 1977. 

Calculated particulate extinction along the flight 
track for one-minute data sets for X= 10.6 urn. 
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Fig. A 10(g). GAMETAG flight data for September 2, 1977. 

Calculated particulate extinction along the flight 
track for five-minute data sets for A= 10.6 ym. 
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Fig. A 10 (i) . GAMETAG flight 
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Table All. Significant times for April 27, 1978. 

Denver, Colorado to Moffett Field, 
California. 


Significant Points 


# 

TIME 

1 

17:00 

2 

17:21 

3 

19:23 

4 

19:24 

5 

19:54 

6 

20:19 

7 

20:25 

8 

20:28 

9 

20:52 

10 

21:14 

11 

22:36 

12 

22:53 


Denver 


Moffett Field 
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Fig. All (a). GAMETAG flight data for April 27, 1978. 

Altitude and location flight track plotted as a 
function of time after takeoff. 
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Fig. AH(d). GAMETAG flight data for April 27, 1978. 

Calculated backscatter coefficient along the flight 
track for five-minute data sets for X = 0.63 ym . 
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Fig. AH (e). GAMETAG flight data for April 27, 1978. 

Calculated ratios for backscatter to extinction for 
five-minute data sets for A= 0.63 urn. 
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Fig. A 11(f). GAMETAG flight data for April 27, 1978. 

Calculated particulate extinction along the flight 
track for one-minute data sets for X = 10.6 pm. 
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DRTE=RPR I L 27. 1978 
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Fig. All (g). GAMETAG flight data for April 27, 1978. 

Calculated particulate extinction along the flight 
track for five-minute data sets for X= 10.6 jam. 
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Fig. All (h). GAMETAG flight data for April 27, 1978. 

Calculated backscatter coefficient along the flight 
track for five-minute data sets for X= 10.6 ym. 
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Fig. All (i). GAMETAG flight data for April 27, 1978. 

Calculated ratios for backscatter to extinction for 
five-minute data sets for A= 10.6 }jm. 
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Table A12. Significant times for May 2, 1978. 

Hilo, Hawaii to Johnston Atoll. 


Significant Points 


# 

TIME 


1 

20:03 

Hilo 

2 

20:19 


3 

20:36 

(Long.) = 20:42 (Alt). 

4 

20:5? 


5 

22:05 


6 

22:28 

(Lat.) = 22:36 (Alt.) 

7 

23:17 


8 

23:48 


9 

00:47 


10 

01:18 


11 

02:21 

Johnston Atoll 
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Fig. A12 (b). GAMETAG flight data for May 2, 1978. 

Calculated particulate extinction along the flight 
track for one-minute data sets for A= 0.63 ym. 
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Fig. A 12 (d). GAMETAG flight data for May 2, 1978. 

Calculated backscatter coefficient along the flight 
track for five-minute data sets for X= 0.63 ym. 
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Fig. A 12 (e) . GAMETAG flight data for May 2, 1978. 

Calculated ratios for backscatter to extinction for 
five-minute data sets for X = 0.63 ym. 
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Fig. A 12(f). GAMETAG flight data for May 2, 1978. 

Calculated particulate extinction along the flight 
track for one-minute data sets for X= 10.6 nm. 
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Fig. A 12(h). GAMETAG flight data for May 2, 1978. 

Calculated backscatter coefficient along the flight 
track for five-minute data sets for X= 10.6 yin. 
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Fig. A 12 (i) . GAMETAG flight data for May 2, 1978. 

Calculated ratios for backscatter to extinction for 
five-minute data sets for A= 10.6 im. 
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Table A 13. Significant times for May 4, 1978. 

Canton Island to Nandi, Fiji Islands. 


Significant Points 


TIME 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 


Canton Island 


20:55 
21:26 
22:37 
23:05 
23:18 

0:11 (Alt.) - 0:22 (Long.) 

0:41 

2:23 

2:45 

2:53 

3:10 

3:26 Nandi 
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Fig. A13 (a). GAMETAG flight data for May 4, 1978. 

Altitude and location flight track plotted as 
function of time after takeoff. 
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Fig. A13 (b). GAMETAG flight data for May 4, 1978. 

Calculated particulate extinction along the flight 
track for one-minute data sets for A= 0.63 yin. 
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Fig. A13 (c). GAMETAG flight data for May 4, 1978. 

Calculated particulate extinction along the flight 
track for five-minute data sets for A= 0.63 ym. 
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***Knol I enberg Data*** 
DRTE=MRY 4. 1978 
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Fig. A13 (d). GAMETAC flight data for May 4, 1978. 

Calculated backscatter coefficient along the flight 
track for five-minute data sets for X= 0.63 ym. 


j 


A125 




RATIO' B K 5 C / 5 I G M A > 


***Kno I I enberg Do+a*** 
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Fig. A13 (e). GAMETAG flight data for May 4, 1978. 

Calculated ratios for backscatter to extinction for 
five-minute data sets for \= 0.63 ym. 
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Fig. M3 (f). CAMETAC flight data for May 4, 1978. 

Calculated particulate extinction along the flight 
track for one-minute data sets for X= 10.6 nm. 
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Fig. A13 (h). GAMETAG flight data for May 4, 1978. 

Calculated backscatter coefficient along the flight 
track for five-minute data sets for X- 10.6 pm. 
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Fig. A13 (i). GAMETAG flight data for May 4, 1978. 

Calculated ratios for backscatter to extinction for 
five-minute data sets for A= 10.6 vin. 
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Table A 14. Significant times for May 10, 1978. 

Christchurch, New Zealand to 

Christchurch, New Zealand. 

$ 

■r 


Significant Points 


TIME 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 


22:08 

22:35 

23:46 

0:01 

0:30 

1:04 

2:40 

3:18 

4:21 

4:28 

4:54 


Christchurch 


Christchurch 
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Fig. A 14 (a). GAMETAG flight data for May 10, 1978. 

Altitude and location flight track plotted as 
function of time after takeoff. 
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Fig. A 14(b). CAMETAG flight data for May 10, 1978. 

Calculated particulate extinction along the flight 
track for one-minute data sets for A = 0.63 ym. 
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Fig. A 14(d). GAMETAC flight data for May 10, 1978. 

Calculated backscatter coefficient along the flight 
track for five-minute data sets for X = 0.63 ym. 
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Fig. A 14(e). GAMETAG flight data for May 10, 1978. 

Calculated ratios for backscatter to extinction for 
five-minute data sets for A= 0.63 ym. 
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Fig. A 14(f). GAMETAG flight data for May 10, 1978. 

Calculated particulate extinction along the flight 
track for one-minute data sets for X-?= 10.6 pm. 
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Fig. A 14 (g) . GAMETAG flight data for May 10, 1978. 

Calculated particulate extinction along the flight 
track for five-minute data sets for A= 10.6 ym. 
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V 

Fig. A 14(h). GAMETAG flight data for May 10, 1978. 

Calculated backscatter coefficient along the flight 
track for five-minute data sets for X= 10.6 ym. 
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***Kno I I enberg Data*** 
DfiTE=MRY 10.1978 



TIME (HOURS AFTER TAKE OFF) 

Fig. A 14 (i) . CAMETAC flight data for May 10, 1978. 

Calculated ratios for backscatter to extinction for 
five-minute data sets for A= 10.6 pn. 




A140 


i LULU 




Table A 15. Significant times for May 11, 1978. 

Christchurch, New Zealand to 
Nandi, Fiji Islands 


Significant Points 


TIME 


1 

2 

3 

4 

5 

6 

7 

8 

9 

★ 

10 


(Alt.) = 00:30 (Long.) 


Computer on out of Christchurch 
22:34 
22:52 
23:50 
00:27 
01:12 
01:33 
03:39 
03:46 
04:03 

Data gap 04:04 - 04:11 
04:20 Nandi 
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Fig. A15 (a). GAMETAG flight data for May 11, 1978. 

Altitude and location flight track plotted as 
function of time after takeoff. 
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Fig. A15 (d). GAMETAC flight data for May 11, 1978. 

Calculated backscatter coefficient along the flight 
track for five-minute data sets for \= 0.63 ym* 
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***Knoi lenberg Date*** 
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TIME (HOURS AFTER TAKE OFF) 

Fig. A15 (e). GAMETAG flight data for May 11, 1978. 

Calculated ratios for backscatter to extinction for 
five-minute data sets for A= 0.63 pn. 
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Fig. Al5(f). CAMETAG flight data for May H. 1978 - 

Calculated particulate extinction along the flight 
track for one-minute data sets for X = 10.6 ^m. 
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Fig. A 1 5 (h) . GAMETAG flight data for May 11, 1978. 

Calculated backscatter coefficient along the flight 
track for five-minute data sets for X= 10.6 ym. 
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Fig. A 1 5 (i) . GAMETAC flight data for May 11, 1978. 

Calculated ratios for backscatter to extinction for 
five-minute data sets for A = 10.6 jm. 
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Table A16. Significant times for May 12, 1978. 

Nartcli, Fiji Islands to Canton Island. 

jr r 


Significant Points 


TIME 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 


21:50 Nandi 

21:56 

22:10 

22:13 

23:11 (Alt.) - 23:36 (Long.) 
23:47 
0:25 
1:56 
2:30 
3:40 

3:45 Canton Island 
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Fig. A16 (a). GAMETAG flight data for May 12 1978 

Altitude and location flight track plotted as a 
function of time after takeoff. 
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Fig. Al6(d). GAMETAG flight data for May 12, 1978. 

Calculated backscatter coefficient along the flight 
track for five-minute data sets for X= 0.63 ym. 
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Fig. Al6(e). CAMETAG flight data for May 12, 1978. 

Calculated ratios for backscatter to extinction for 
five-minute data sets for A= 0.63 pm. 
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***Kno I 1 enberg Data*** 
DRT E =MRY 12. 1978 



Fig. A 16 (h) . GAMETAG flight data for May 12, 1978. 

Calculated backscatter coefficient along the flight 
track for five-minute data sets for A= 10.6 ym. 
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Fig. A 16 (i) . CAMETAC flight data for May 12, 1978. 

Calculated ratios for backscatter to extinction for 
five-minute data sets for X= 10.6 im. 


i 


A160 




Table A 17. Significant times for May 14, 1978. 

Johnston Atoll to Hilo, Hawaii. 


Significant Points 


# TIME 


1 

2:00 

Johnston Atoll 

2 

2:32 


3 

3:27 


4 

3:37 


5 

3:44 


6 

4:05 (Alt.) - 4:03 (Long.) 

7 

4:45 


8 

4:54 


9 

5:36 


10 

6:32 


11 

6:38 


12 

7:01 


13 

7:35 


14 

7:45 


15 

7:59 

Hilo 
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Fig. A 17 (a) . GAMETAG flight data for May 14, 1978. 

Altitude and location flight track plotted as a 
function of time after takeoff. 



*1X3 





EXT. 





BKSC . X l.L + 7 (1/M . 5TR) 


ORIGINAL PAGE IS 
OF POOR QUALITY 


* * * K n o i lenberg Dat'd*** 

DRTE=MRY 14.1978 

I 


TIME (HOURS Ri" T ER TAKE Oi-F) 


Fig. A 17 (d) . GAMETAG flight data for May 14, 1978. 

Calculated backscatter coefficient along the flight 
track for five-minute data sets for X= 0.63 yin. 
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Fig. Al7(e). CAMETAG flight data for May 14, 1978. 

Calculated ratios for backscatter to extinction for 
five-minute data sets for X= 0.63 ym. 
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Fig. A 17 (f J . GAMETAG flight data for May 14, 1978. 

Calculated particulate extinction along the flight 
track for one-minute data sets for X= 10.6 pm. 
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Fig. A 17(h). CAMETAG flight data for May 14, 1978. 

Calculated backscatter coefficient along the flight 
track for five-minute data sets for X = 10.6 yin. 
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Fig. A 17 (i) . GAMETAG flight data for May 14, 1978. 

Calculated ratios for backscatter to extinction for 
five-minute data sets for A= 10.6 nn. 
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Table A18. Significant tiems for May 18, 1978. 
Moffett Field, California to 
Denver, Colorado. 


Significant Points 


TIME 


1 

2 

3 

4 

5 


17:56 

18:21 

20:23 

20:40 

20:57 


Moffett Field 


Denver 
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Fig. A 18 (a) . GAMETAG flight data for May 18, 1978. 

Altitude and location flight track plotted as 
function of time after takeoff. 
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TIME (HOURS AFTER TAKE OFF) 


Fig. A 18 (b) . GAMETAG flight data for May 18, 1978. 

Calculated particulate extinction along the flight 
track for one-minute data sets for X= 0.63 ym. 
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Fig. A 18(d). GAMETAG flight data for May 18, 1978. 

Calculated backscatter coefficient along the flight 
track for five-minute data sets for A= 0.63 ym. 
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Fig. A 18(e). GAMETAG flight data for May 18, 1978. 

Calculated ratios for backscatter to extinction for 
five-minute data sets for X= 0.63 vtm. 
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Fig. A 18(f). CAMETAG flight data for May 18, 1978. 

Calculated particulate extinction along the flight 
track for one-minute data sets for X= 10.6 pm. 
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Fig. A 18(h). GAMETAG flight data for May 18, 1978. 

Calculated backscatter coefficient along the flight 
track for five-minute data sets for X= 10.6 ym. 
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Fig. A 18 (i) . GAMETAG flight data for May 18, 1978. 

Calculated ratios for backscatter to extinction for 
five-minute data sets for A= 10.6 ym. 
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Table A19. Significant tiems for May 27, 1978. 

Denver, Colorado to Great Falls, 
Montana,. . 


Significant Points 
# TIME 

Denver 


Great Falls 


1 

2 

3 

4 

5 

6 


22:04 

22:44 

23:44 

0:00 

0:13 

0:21 
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EXT. X l.E+7 (1/M) 


***Knol 1 enberg Dq+q*** 
DATE = MA Y 27. 1978 



TIME (HOURS AFTER TAKE OFF) 

Fig. A 19(c). GAMETAG flight data for May 27, 1978. 

Calculated particulate extinction along the flight 
track for five-minute data sets for A= 0.63 ym. 
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DflTE=MflY 27. 1978 
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Fig. A 19(d). GAMETAG flight data for May 27, 1978. 

Calculated backscatter coefficient along the flight 
track for five-minute data sets for X = 0.63 ym. 


4 


A185 



RAT'D (BKSC/SIGMfl) 
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Fig. A 19(e). GAMETAG flight data for May 27, 1978. 

Calculated ratios for backscatter to extinction for 
five-minute data sets for X= 0.63 pm. 
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EXT. X l.E+7 (1/M) 
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***Kno I lenberg Goto*** 
DATE=MRY 27.1978 



Fig. A 19(f). GAMETAG flight data for May 27, 1978. 

Calculated particulate extinction along the flight 
track for one-minute data sets for X= 10.6 tim. 
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***Knol I enbsrg Do+o*** 
DRTE=MflY 27.1978 



TIME (HOURS RFTER TRKE OFF) 


Fig. A19 (h). CAMETAG flight data for May 27, 1978. 

Calculated backscatter coefficient along the flight 
track for five-minute data sets for A= 10.6 pm. 
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RATION* 5C/SIGMA) 


***Knol I enberg Data*** 
DflTE=MflY 27.1978 



Fig. A 19 (i) . CAMETAG flight data for May 27, 1978. 

Calculated ratios for backscatter to extinction for 
five-minute data sets for X= 10.6 
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Table A20. Significant times for May 28, 1978. 

Great Falls, Montana to White Horse, 
Canada . 


Significant Points 


# TIME 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 


16:09 Great Falls 

16:31 

17:05 

17:09 

17:22 

17:46 

18:45 

19:11 

19:52 

21:03 

21:22 White Horse 
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Fig. A20 (a). GAMETAG flight data for May 28, 1978. 

Altitude and location flight track plotted as 
function of time after takeoff. 





EXT. X 1 . E + 7 <1/M> 


fc+^Knol I enberg Data*** 
DRTE=MRY 28. 1978 



Fig. A20 (b). GAMETAG flight data for May 28, 1978. 

Calculated particulate extinction along the flight 
track for one-minute data sets for \= 0.63 pm. 
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EXT. X l.E + 7 ( 1 / M ) 


***Kno i I enberg Do+o*** 
DRTE=MRY 28.1978 



TIME (HOURS RFTER TRKE OFF) 


Fig. A20 (c). GAMETAG flight data for May 28, 1978. 

Calculated particulate extinction along the flight 
track for five-minute data sets for X= 0.63 ym. 
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***Knoi I enberg Da+o*** 
DRTE=MRY 28.1978 



Fig. A20 (d). GAMETAG flight data for May 28, 1978. 

Calculated backscatter coefficient along the flight 
track for five-minute data sets for X= 0.63 ym. 
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* * * K n o I lenberg Data*** 
DflTE=MflY 28. 1978 



TIME (HOURS RFTER TRKE OFF) 

Fig. A20 (e). GAMETAG flight data for May 28, 1978. 

Calculated ratios for backscatter to extinction for 
five-minute data sets for A= 0.63 urn. 
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EXT. 


***Knol ! enberg Goto*** 
DATE =MAY 28.1978 



TIME (HOURS AFTER TAKE OFF) 


Fig. A20 (f). GAMETAG flight data for May 28, 1978. 

Calculated particulate extinction along the flight 
track for one-minute data sets for X= 10.6 nm. 
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EXT. X t.E+7 (1/M) 


***Knoi I enberg Data*** 
QRTE=MflY 28. 1978 



TIME (HOURS AFTER TAKE OFF) 


Fig. A20 (g). GAMETAG flight data for May 28, 1978. 

Calculated particulate extinction along the flight 
track for five-minute data sets for A= 10.6 yin. 
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***Knol I enberg Data*** 
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Fig. A20 (h). GAMETAG flight data for May 28, 1978. 

Calculated backscatter coefficient along the flight 
track for five-minute data sets for A= 10.6 um. 
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Fig. A20 (i). GAMETAG flight data for May 28, 1978. 

Calculated ratios for backscatter to extinction for 
five-minute data sets for A= 10.6 jin. 
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Table A21. 


Significant tiems for May 31, 1978. 
White Horse, Canada to Great Falls 
Montana . 


* ** ,* 


Significant Points 


# 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 


TIME 

18:14 White Horse, Canada 
18:41 
18:53 
19:21 
19:36 
20:09 
20:51 
21:04 
22:21 

22:43 Great Falls 
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Fig. A21 (a). CAMETAG flight data for May 31, 1978. 

Altitude and location flight track plotted as 
function of time after takeoff. 
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Fig. A21 (b). GAMETAG flight data for May 31, 1978. 

Calculated particulate extinction along the flight 
track for one-minute data sets for X= 0.63 pm. 
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Fig. A21 (c). GAMETAG flight data for May 31, 1978. 

Calculated particulate extinction along the flight 
track for five-minute data sets for X= 0.63 ym. 
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Fig. A21 (e). GAMETAG flight data for May 31, 1978. 

Calculated ratios for backscatter to extinction for 
five-minute data sets for X= 0.63 yin. 
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Fig. A 21 (g) . GAMETAG flight data for May 31, 1978. 

Calculated particulate extinction along the flight 
track for five-minute data sets for A= 10.6 ym. 
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Fig. A 21 (h) . GAMETAG flight data for May 31, 1978. 

Calculated backscatter coefficient along the flight 
track for five-minute data sets for X = 10.6 ym. 
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***Knol I enberg Dota*** 
DRTE=MRY 31. 1978 



Fig. A 21 (i) . GAMETAG flight data for May 31, 1978. 

Calculated ratios for backscatter to extinction for 
five-minute data sets for A= 10.6 im. 
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Table A22. Significant times for June 1, 1978. 

Great Falls, Montana to Denver, 
Colorado. 


Significant Points 
# TIME 

Great Falls 
Denver 


1 

2 

3 

4 

5 

6 


15:50 

16:05 

16:37 

16:44 

17:36 

17:45 
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Fig. A22 (a). GAMETAG flight data for June 1, 1978. 

Altitude and location flight track plotted as a 
function of time after takeoff. 










TIME (HOURS RFTER TRKE OFF) 


Fig. A22 (b). GAMETAG flight data for June 1, 1978. 

Calculated particulate extinction along the flight 
track for one-minute data sets for A= 0.63 ym* 
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Fig. A22 (d). GAMETAG flight data for June 1, 1978. 

Calculated backscatter coefficient along the flight 
track for five-minute data sets for X - 0.63 pm. 
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***Kno I I enberg Goto*** 
DflTE=JUNE 1.1978 



Fig. A22 (e). GAMETAG flight data for June 1, 1978. 

Calculated ratios for backscatter to extinction for 
five-minute data sets for X= 0.63 ym. 
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Fig. A22 (h). GAMETAC flight data for June 1, 1978. 

Calculated back scatter coefficient along the flight 
track for five-minute data sets for X = 10.6 ym. 
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Table B1. Combined Size Ranges for GAMETAG 



Optical Particle Counter Data 


Instrument 

Range # 

1977 

Radius Interval (pm) 

Mean 

Radius (pm) 

ASAS 

1 

0.126 - 0.170 

0.146 

ASAS 

2 

0.170 - 0.274 

0.216 

ASAS 

3 

0.274 - 0.385 

0.325 

ASAS 

4 

0.385 - 0.500 

0.439 

ASAS 

5 

0.500 - 0.625 

0.521 

ASAS 

6 

0.625 - 0.750 

0.685 

FSSP 

7 

0.750 - 1.104 

0.91 

FSSP 

8 

1.10 - 1.73 

1.38 

FSSP 

9 

1.73 - 2.34 

2.01 




1978 



ASAS - 

Same as 1977 




FSSP 

7 

0.25 

- 0.75 

0.43 

FSSP 

8 

0.75 

- 1.73 

1.14 

FSSP 

9 

1.73 

- 2.34 

2.01 
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Table B2. Size Data (dn - cts/cm 3 ) for 1977 Flights (continued) 
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Table B2. Size Data (dn - cts/cm ) for 1977 Flights (continued) 
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Table B2. Size Data (dn - cts/cm ) for 1977 Flights (continued) 



CO 

m 

c 

c 

C 



1 

w 

u. 

UJ 

o 

C' 

ir 

c a 

c a 

P« | 

c 1 

C 1 

is a 

• i 

• a 

• t 

a 

a 


a 

i 


a 

a 


a 

a 


ri 

r- i 

cr i 

c a 

c a 

C 1 

a a 

i a 


U. 1 

UJ 1 

kU 1 

-X | 

c a 

Pm | 

-r a 

pj a 

P* I 

— i 

04 | 

O' 1 

• a 

•1 

* a 

a 

i 

l 

a 

a 

l 


— a 

t 

** i 

04 ( 

C 1 

c a 

c a 

a a 

a i 


w a 

kU 1 

um a 

c a 

p* a 

p i 

p*a 

^ a 

C 1 

IT 1 

u 1 

04 | 

• 1 
1 
1 

• a 
a 
a 

1 | 

1 

1 

— 1 

a 

— i 

c a 

C t 

a i 

c I 

l 1 
kU 1 

• i 
Ud 1 

uj a 

c 1 

r 1 

od t 

S' 1 

p- 1 

-r » 

m 1 

— 1 

PM | 

• a 

• 1 

• | 


I 

I 

I 

I 

C I 
I I 
U. I 
N I 
^ I 

pj | 

• i 
I 

I 

l 

t 


I 

l 

I 

* 

o i 
l t 
U I 
-r t 
IT | 
f 

• | 
I 
« 
I 
I 


I 

I 

I 

t 

C I 

O I 
I 

u. I 
r- I 
%r » 
is; i 
• i 
i 
i 
I 


— 1 

_ | 

c t 

c* 

c 1 

c a 

c a 

c 

1 1 

a a 

i 


k~ 1 

u> a 

m 1 

u. 

e* 1 

a 

«»• 1 

r- 

p* a 

m J 

n i 

c 

pm i 

PJ | 

pj » 

a- 

• a 

O I 

• a 

• 


• i 

i 

i 

i 

I 

c i 

l 

U. I 
C I 
IT I 

«r I 

• i 

» 

1 

I 

i 

«-4 | 

n i 

I 

u I 
C » 
c: i 

IT * 

• I 
I 
I 
I 
I 

C I 
e i 
i 

UJ I 

r* i 
or i 

• i 

I 

i 

i 

i 

H | 

C I 

I 

U t 
C I 
C I 

o*4 | 

• | 
I 
I 
I 
I 


G t 

c a 

m I 

n a 

o a 

o a 

O 1 

o a 

a 

i 

1 


UJ I 

UJ I 

uj a 


*»• a 

c a 

m a 

04 | 

u> a 

m | 

m a 

P J * 

C 1 PCM 1 

cr-'M i 

OCC 1 

• a 

mm • a 

mm • a 

mm • a 


w 

© 


P4 • 

• I 

I 

I 

I 

a 


a 

i 

i 

i 

o i 
i 

v*j a 

N | 
IS I 
If* I 
© | 
I 


I 

c 1 
c » 
a 

uj I 

r f 
m l 
tr j 

• i 

t 

l 

I 

l 

C I 
C I 
I 

u I 

—• t 

r* i 

r- i 

• i 
i 
» 
i 
i 

c i 

© i 

i 

i-» i 

p* t 

• t 


O 1 

m | 

C 1 

— i 

o4 | 

e a 

C l 

c a 

r* i 

a i 

a 

1 1 

a 

i 

a 

uj a 

UJ l 

UJ I 

UJ 1 

ui a 

c » 

pj 1 

o4| 

bX 1 

p i 

m i 

^ 1 

'‘J 1 

«p i 

Pi a 

n 1 

p- a 

r a 

pj a 

pj a 

• 1 

• a 
i 

• a 

• i 

• » 


a 

a 

a 


c; 

c 


PJ | 
— I 


I 

I 

I 

I 

04 | 

C I 

a 

u. i 

C I 

>h I 

t 

• i 
t 

a 

i 

» 

© I 
O t 
I 

L. I 
r*- | 
-* I 
c* i 

• i 

l 

i 

i 

» 

r* i 

» 

u. t 
w* » 

•4 | 

• | 

I 

t 

a 

t 

m 1 

c i 

i 

k-J t 

S. I 

C I 

-4 | 
• | 
I 



u. I 
p* I 
*0 I 
CM I 
• I 

•a | 

o a 

VJ I 

r* i 

*>■ t 

rc cicr «r i 

• i • i 

•• •• i •••• a 
so i sp* » 
t.*w a <h»-. i 
•••• a a 

0 - i oo ^ a 
*mpj © a isr'J o a 
n n a ■ n a 

o I -Jl 

1 - ' I uuC I 

-r i *.*. o a 

— — m l ■ — ^ a 

k-b- • I M- • 8 


| 

• | 
a 
a 
a 
I 


u- • 

c* » 

• a 


V? 

C I 
C I 
• • 
a 
a 
a 


c 

c 


u 

c 


m 1 

m a 

P* | 

40© 1 

pj a 

m 


m 


mm | 

P. | 

*4 


m 

1 

IV 


O 1 

U 1 

C I 

C 1 

c: a 

c 


c 


Q 9 

c a 

c 


C 

1 

c 


1 

1 

a 

a • 

a i 

l 


l 


1 1 

I 




1 

1 


U. 1 

W 1 

w f 

u 1 

y 1 

Ui 


UJ 


W 1 

u 1 

k. 


Ui 

a 

UJ 


c a 

»r a 

r a 

p» t 

r a 

« 


-X 


t* 1 

C 1 

•X 


e 

a 

c 


G 1 

r a 

o4 | 

m i 

pj a 

o 


p- 


m a 

bp 1 

«c 


*c 

i 

p- 


*r i 

-<• a 

P* 1 

<T l 

r. i 

r 


m 


04 | 

04 | 

o4 



• 

fiO 


• a 

• a 

• | 

• | 

• a 

• 


• 


• | 

• | 

• 


• 

a 

• 



a 

a 

a 

o4 | 

O I 
a 

uj a 

N | 

«4 | 


r i 
m t 

m a 

PJ 

I* 

ki 


a 

a 

a 

a 

c* a 

a 

u a 
c i 
tr. | 

• i 

a 

a 

a 

a 

i a 
it a 

• i 

a 

a 

a 

a 

«— i 

r* i 
a i 

i 

c* » 

-r i 

• a 

i 

a 

a 

a 

r a 
r a 
l a 

u. a 

•r i 
r. i 
C* I 

• I 

a 

a 

a 

a 

o a 
w a 

Pi I 
P* | 

lb | 

• I 

a 

a 

a 

a 

c a 
© a 
a 

cj a 


e a 
c • 


*r l 

i 


r 

c 


a 

a 

a 

i 

c a 
a a 

U I 

c a 
? a 
*r i 

• i 
a 
i 
a 
I 

c a 
a a 

UJ I 

r- a 
m t 

• a 
a 
a 
a 
a 

a 
a 

pj a 


c 

a 


V m 


OVJ 

1 oo 

a oo 

| OO 


« z? 

» rr 

a «r 

■no 

1 MM 

a 

a 

►-a 

I WQ 

1 mC 

1 me 


: *a a a. ** « a at 
+ I <u | 

:: a *- z t *- 
c i !/•» OIlK 



. , . 'ICCP 

I r*n • i mm • 

c* a •• •• a •• •• 

• i „*■*» t> a po 
tr I -r a p 
•j i •• •• a •• •• 

• | Oo- Pd a 

- | <bj«* i — * | G 

y.i m a *a u 

a w i u 

• a uU'iT a kju.iv 
i a i j.> — 

l i a *— a 
j a mm • a m*- 
— a a 

m a oo a jo 
a rr' a rr* 

*t a — - a 
»- a me a ► n 


cc-' 

an • 


cap* 

pg*r 


pjpj© 
«a h 




OO 


a 

a 

a 

a 

a h 

a o 

a 

a uj 

a ip 

a r> 

a r CbO 

I nn • 

I •• •• 

I ?« 

I o 

I Hfwfg 

a ivjimo 

| N N 

a u 

a uuip 
irvN 
a 

I mm 

I 

a oo 
t rr 



4r 

C 

a 

Uj 

r 

Pm 


C 


c 

a 


o 

a 


i 

Pi | 
• | 
I 
I 
I 


© 

a 

u a 
p* a 
— a 
<£ 


a 

a 

a 

a 

c a 
a a 

u. a 

c a 
r a 

• I 

I 

a 

a 

a 

c a 
c i 
i 

Uj a 
c a 
© a 
g a 

• a 
* 
a 
a 
i 

pj a 
c a 
a a 

u. I 
si a 
pj 


a 

a 

a 

a 

— a 
c 


a 

i 

a 

a 

h a 
o a 
a 

u* a 
■«* i 

m4 | 
PJ I 
• | 
a 
a 
a 
a 

c i 
o a 
i 

w a 


bf a 

m 1 

<M | 

C* 1 

m a 

m i 

r* a 

m | 

• i 

• a 

• | 

• 1 

• i 

* a 

• a 

• a 

a 

a 

• 

1 

a 

• 

a 

i 

a 

1 

a 

» 

t 

I 

1 

■ 

* 

a 

a 

■ 

a 

i 

a 

1 

1 

0-4 | 

1 

—4 | 

e a 

i 

c a 

* 

C 1 

— i 

C 1 

c a 

C 1 

c a 
■ • 

C 1 

i 

o a 
• 

C 1 
1 

O 1 

a a 

kJ 1 

UJ I 

UJ 1 

i • 
UJ I 

UJ • 

UJ 1 

» 

Uj 1 

Uo 1 

m a 

O 1 

C | 

bC a 

pj » 

c a 

mo a 

c* a 

*0 a 

04 | 

p*b 1 

•4 | 

O 1 

n a 

C | 

p- a 

bX 1 

tb a 

r a 

pj a 

n i 

m a 

Pi 1 

bC 1 

a 

• • 
a 

• 1 
a 

• a 
a 

• a 
a 

• i 

a 

* a 
• 

• a 

a 

a 

a 

a 

a 

a 

a 

a 

i 

a 

a 

» 

a 

a 

» 

a 

a 

a 

a 

a 

o a 

c t 

c a 

o4 | 

o a 

o a 

C 1 

c a 

c a 

o I 

O 1 

C 1 

o I 

O 1 

o a 

o » 

a 

a 

a 

1 a 

a 

i 

a 

i 

uj a 

UJ 1 

UJ a 

UJ 1 

UJ 1 

uj a 

Uj | 

UJ 1 

K | 

m a 

C_- | 

O | 

o a 

o4 | 

•4 | 

Pw | 

•r a 

c- a 

PJ 1 

— 1 1 

o t 

o a 

tr j 

pj i 

Pd i 

Po a 

«4 | 

p 1 

C 1 

or* | 

P 1 

pj a 

• a 

• a 

• | 

• a 

• a 

• | 

• i 

• i 


i 
a 
a 
a 

-4 | 

O ■ 


ar 

T9 


a c a 

a a 
ui a uj i 

«? i -• * 

ci i- i 

ec«r i cc-* i ccpj a 

men • | nn • a r*ry « | mm • 

•• •• | •• •• | •• •• | •• •• 

o? | — <0 I PJ*n | mrsi 

mru I rntP I pm | mm 

•• •• I | •• •• | •• •• 

psivo | Pjiv.m | fsmm | nnn 
PJPJO a PjPjO I PJPJO I pJpjO 
■ n 100 | 0 0 I 0 0 

j a uj a j i uj 

uwP I UjujO | ujUi«t I uu-O 



x scg t a 

| m t — - .. 

I • l nn 



kP O I kP o 

a 

a 

a 


t oo 

I OO 

t oo 


a err 

a if e 

1 M44 

| 0404 

1 04-4 

a me 

1 me 

1 mo 

1 OCdC 0 

1 JCj: a 

a «2r 

a <u. 

I «UJ 

1 <Ui 

1 m ^ 

i m rr 

1 m 

1 O 

a ui ca 

a kP i 


a 

a 


a 

a 


a 

a 


a 
a 
a 

i 

— a 
o a 
a a 

UJ i 

p a 

r a 
IT I 

• a 
a 
i 

a 

l 

m | 
O I 
I I 

UJ I 

<" I 
| 

— a 

• a 

» 

i 

a 

— i 

c a 
t a 
-j a 
<-? » 
a 


a 

a 

a 

o i 

o a 
» 

UJ I 

N I 
o a 
cc-c a 

nn • i 

• 4 •• | 

a 
» 
a 

r* 1 | 

pjpjO I 





B7 



Table B2. Size Data (dn - cts/cnr ) for T977 Flights (continued) 


© 


O 

U 

o 

PVI 




00 

U 

© 

o 

©CO 

-0 

c 


c 

c 

o 

C 

© 

© 


© 

© 

© 

© 

o © c 

c 






• 

1 

1 







1 

l*J 


tti 

UJ 

u. 

UJ 

UJ 

UJ 


u 

Ua 

U 

Ut 

u u u. 

UJ 

o 


c* 

u 

c 

© 

© 

m 



U 

© 

© 

© © CJ 

<£ 

o 


c e 

O 1 

O 1 

r t 

© 1 

c I 


© 9 

© # 

© 1 

O 1 

Cl C| Cl 

m 

c 


c I 

O 9 

C 9 

m 1 

«V 1 

ar 9 


ru 9 

© • 

© 1 

c » 

© 1 ©1 © 9 

9 *- 

0 


• i 

# | 

0 | 

0 | 

0 1 

0 i 



0 | 

0 | 

0 I 

0 I 0 1 0 | 

1 1 1 

9 1 1 

0 

© 


r» f 

r | 

fSi | 

C 9 

P 9 

p i 


•0 1 

fVl 1 

fo | 

r i 

9 1 9 

1 1 1 

1 ©| M| 

© 

o 


C 1 

C 1 

O 1 

C 9 

C 1 

O 9 


© 1 

© 1 

© 1 

C 1 

C 9 © 1 © 1 

© 

1 


1 i 

1 1 

1 1 






9 1 

9 1 

1 1 

9 9 9 1 9 1 


Uj 


UJ I 

u- 9 

Uj 9 

u J 9 

Uj 9 

uj 9 


UJ 1 

i.j • 

U | 

W 1 

Uj 1 Uj | w 1 

U 

© 


C 1 

C 1 

<M 9 

© 1 

© 1 

© 1 



© 1 

C- 1 

if 1 

m 1 m | mi 


it 


•C 1 

CM 1 

C 1 

Pj | 

*C | 

© 1 


N 1 

© 1 

o 9 

© i 

m | 00 1 p i 

© 

IS* 


O' 1 

© 1 

N 9 

n i 

p 9 

p 9 


® l 

© 9 

m i 

r* i 

© 1 n i P| 

© 

• 


• | 

* 1 

o 9 

0 9 

0 9 

0 9 






0 | • j 0 j 

f I I 

1 1 1 

0 

•4 


00 | 

H | 


cr 9 




1 


© 1 

„i 

1 9 1 

1 t 1 

— 1 © 1 © 1 

© 

o 


C 1 

O 1 

© 1 

© 9 

O 1 

C 1 


O 9 

O 9 

O 9 

O | 

© 1 © 1 O 1 

0 

1 


I f 

9 1 

9 9 






1 9 

1 9 

9 1 

19 II 1 


w 



u. 1 

U* | 

uj i 

u | 

Uj 9 



Uj 9 

U. | 

© 1 

Uj 1 UJ | Ul 

u 

m 


CC 1 

p- 1 

CM 1 

*C t 

O’ 9 

C 1 


h> | 

© 9 

C 1 

O 1 

i © i mi 

O 

rc 


m 1 

m » 

p 1 

© t 

C I 

© 9 


© 1 

© 1 

© 9 

4-0 1 

n- 1 © i © i 

m 

«T 


P- 1 

h i 

P- 1 

IT I 

rv. | 

© 1 


N 1 

© 9 

1 

ar i 

© I © 9 0-9 

© 

• 


• | 

# i 

0 | 

0 | 

0 | 

0 | 






0 9 • | 0 | 

9 1 1 

1 1 1 

0 



CM 1 

r 4 | 

C 1 

*0 | 

© 1 

C | 


M 4 | 


© 1 

. J 

9 9 9 

1 1 1 

© 1 M p 1 

© 

O 

1 

u 


o 1 

C 1 
t 1 

t_ } 

O 1 

© 1 

O 1 

© * 


© 1 

U 1 

O 1 

© » 

© 1 © 9 © 1 

O 


UJ 1 

u i 

Uj 1 

U 1 

U 1 


Uj 1 

U 1 

uj i 

UJ 9 

Ui » Ul ul 

UJ 



e i 

© i 

© » 

ro t 

p* i 

r* I 



m 9 

•a. 9 


09 © i mi 

0- 



O 1 

r\j i 

C 1 

C 1 

«£ 1 

X t 


© 1 

m i 

© 1 

Oj | 

mi © i xi 

0 j 

© 


n | 

h t 

© 1 

O' 1 

O' 1 

m i 


O* 1 

00 1 

<Nl 9 

00 1 

p i pi 0-i 

CL 

• 


• i 

• i 

0 1 

0 1 

0 i 

0 i 


0 f 

0 f 

0 | 

0 f 

0 1 0 f 0 1 

1 1 1 

9 9 9 

0 



aJ 1 

PS. 1 

«M 1 

© 1 

C 1 

© 1 


! 

I 

n. t 

J 

1 9 9 

1 9 9 

c 

O 


© t 

O 1 

© 1 

© 1 

© 1 

© I 


© 1 

C 1 

p » 

© 9 

© 1 O « ©1 

© 

1 


1 1 

1 f 

) » 





1 • 

9 1 

9 9 

1 1 

II 19 II 


Uj 


Uj 1 

u 1 

U. 1 

UJ 1 

Uj 1 

Ui t 


UJ 1 

UJ 1 

© 1 

© 1 

Uj 1 U 1 UJ 1 

UJ 

fp 


«o 1 

1 

C 1 

r. i 

•r | 

«M | 


OJ 1 

C 1 

O | 

*T I 

© i tri mi 

O' 




oj 9 

«M 1 

© i 

M | 

© 1 


© 9 

© 1 

O- 1 

m 9 

Ol 1^ t mi 

© 



f\j 1 

© 1 

© t 

p 1 

«r i 

© 1 


aT 9 

\ 1 

© 9 

© 9 

mi mi 0-1 


« 


• ) 

0 | 

0 t 

• 9 

0 9 

0 t 



0 1 



• | 0 f • | 

1 1 1 

1 1 t 

• 




0— | 

| 

C 1 

© 1 

© 1 


0- 1 

aM 1 

© 9 

© 1 

1 9 1 

1 1 1 

c 

c 


O 1 

C 1 

O 1 

C 1 

O 1 

© 9 


© 1 

O 9 

O | 

© 9 

© 9 © 9 O 1 

O 

UJ 


UJ t 

u. t 

Ul 1 

UJ | 

u 1 

U. 1 


UJ 1 

Ul 1 

UJ 9 

Ul 1 

UJ 1 <JJ 1 U 1 

UJ 

m 




© 1 

N | 

< I 

© t 


© 1 

© 1 

<C 1 

© 9 

0.'l CM | oi 

© 

05 


»M 1 

H | 

O 1 

•H | 

rv 1 



© 1 

<C 1 

© 1 

© 1 

»l ro | 0-| 

ee 

•“4 


© 1 

no t 

p 1 

^0 t 

© I 

© 1 


ar i 

m 9 

M I 

r | 

m i mi rv i 

X 

• 


• | 

0 | 

0 t 

0 | 

0 | 

0 t 


0 i 

0 | 

0 | 

0 9 

0 | 0 | • | 
• 1 9 

1 l I 

0 



0-* | 

C 1 

•0 j 

© 9 

© 1 

p 1 


© l 

© | 

© 1 

00 1 

9 i 1 

1 1 9 

©1 © 9 © 1 

p 

O 


C I 

O 1 

O 1 

O 1 

O 1 

a i 


O f 

© 1 

O 1 

© 1 

© 1 Ol © 1 

0 

1 


1 1 


1 1 







9 9 

9 1 

1 1 9 1 


UJ 


U-J 1 

UJ 9 

UJ 1 

UJ t 

UJ 1 

UJ 1 


mJ 9 

UJ 1 

UJ t 

UJ 1 

IP 1 UU 1 UJ 1 

UJ 



O 1 

m i 

© 1 

O 1 


rw | 


P 9 

© 9 


*0 l 

mi mi c » 

*0 

m 


•4 | 

*■0 | 

© i 

m i 


K | 



© 1 

© 1 

m i 

mi mi © i 

CO 

f- 


P» ! 

*0 1 

© 1 

© | 

O i 

P 1 


fi -ri 

© 9 

© 1 

IT | 

CP 1 P. | P» | 

© 

• 


• | 

* 9 

0 9 

0 | 

0 I 

0 | 


9 9 • » 

9 9 1 

1 1 9 

0 | 



0 | 0 | 0 | 
9 9 1 

9 9 1 

0 

o 


O t 

O 9 

C • 

O | 


mJ | 


1 1 1 

1 1 9 

9 9 00 1 

O 9 

O f 

O 9 

1 9 1 

9 1 ! 

© 1 © 1 © 1 

© 

o 


O 1 

O 1 

O 9 

O 9 

O 9 

O 1 


9 9 O t 

O 1 

O 1 

© 9 

© 1 Cl O 1 

0 

w 


U) s 

UJ 9 

UJ • 

UJ 1 

UJ | 

UJ | 


19 UJ 1 

Ul 9 

UJ 1 

Ul | 

Ul 1 ui | Ul 1 

UJ 

IP 


m i 

O t 

n i 

pj t 

^ 1 

cc t 


• 9 fill 

U? 1 

m i 

p • 

mi mi mi 

X 

ru 


—t 1 

O 1 

© 9 

o | 

-C 1 

p I 


19 © | 

m i 

ro i 

m 9 

Ol ro 1 © I 

r\j 

0 ©fNJ 


cor) 9 

CCCI 

CCn 1 

©c© i 

©CP' 1 

©C© 1 

© 

1 i rcfSj i 

©C©> 1 

OC© 1 

cc© 1 

OO© 1 o©m | ©©00 | 

c»©n* 

nw • 


© ro ■ i 

np 0 1 

OP 0 | 

©© 0 1 

© © • 1 

OO t | 

© 





pp • | na 0 | nn 0 1 

©© 0 

aa as 





aa aa 1 

aa aa | 

aa aa | 

aa 

9 | aa aa | 

a* aa | 

aa aa | 

aa aa | 

aa aa | aa aa | aa a* | 

aa aa 

om 


r- o | 

aois. | 

O CO | 

O? 1 

HO | 

OJ -0 1 

© 

1 © I mm | 

*om 9 

^«o 1 

cjr- i 

©V I O© | ©0 1 

r\j~« 

0—#P 


© m { 

in © l 

—j© | 

mm i 

M 1 

| 

«r 

9 0 1 m • 

ISj | 

r j m i 

«r 1 

ng | 00m | m© 1 

ap© 

•• •• 




•• •• i 

aa aa | 

aa aa | 

aa aa 9 

<a 

1 00 I aa aa | 

aa aa | 

aa aa | 

aa aa | 

•a aa | aa aa | aa •« j 

•« aa 

^'I , 0 


■ 4 ’*-! 1 


mine 9 

ihlTO 1 

C« 0 © 5 

©<CfM 1 

o 

« • 9 ©COj | 

©jOO 9 

ooo f 

0 ©c 9 

©©© 1 PPP | PMO f 

ojrg-0 

NNO 


OjPjO 1 

njiMO 1 

no no© | 

OJJNlO I 

INjOJO I 

iNjnJO 1 

rsJ 

1 Z 9 PMO 1 

rsjfMO | 

fSJ«NlO 9 

MNO 9 

NNO | ngnjo 9 CNjrg© 1 

ojrsjO 

M It 


ft N I 

ft « i 

n n 9 

0 ft 1 

N 0 | 

N M | 

M 

lafi h t 9 

ft ft 9 

ft ft | 

ft ft 1 

ft ft 1 ft II | ft ft | 

ft ft 



LJ 1 


UJ 9 

Ul 1 




19 UJ | 

UJ t 

ui 9 

UJ 1 

UJ 1 Ul 9 Ul | 

Ul 

Uku<f 


LJUjO t 


UIUCT 1 

uiu:© | 

lJUP | 

wwm 9 


I ft i uju;m 9 

uut*' t 

UIUJ 00 1 

uiujm 9 

UIUI© | UU'M | UUO 1 

UUJ'O 

».i.X 


LX.O 1 

i 

1 JuC* I 

ISO 1 
>-*— 0 f 

f a;® i 

r;ro | 

a: 

9 1 fc_X N 1 

1 U 9 — — © 1 
* -J 1 ©H- 0 I 

va.®i 

X- Ju© I 

L UO 1 

■Kb.mi r 2 - © 1 a- xx 1 

i_*_m 

P-P* e 


p- p- 0 | 

fc-u • » 

HU 0 1 

uu% 5 

PU 0 f 

w— 

©© 0 i 

©© 0 i 

up 0 i 

•-P • i pp 0 i hp • i 
9 t 1 

P-H- 0 


1 — 1 9 

oo 


OO | 

019 1 

190 1 

OO 9 

OO | 

OO 1 

o 

9 U. 1 OO | 

OO | 

OO | 

OO 9 

00 1 00 9 00 1 

OO 

7 2 1 


•» ■* i 

ie-sr » 

7 r 9 
1 

Z 2 T 1 

MM | 

z ;* » 

7 7 9 

.r 

M 

9 1 if Z 1 

| mx I . 1 

."!? » 
PM | 

;T 7 9 

MM 1 

77 9 

PM | 

72 9 77 9 7? 9 

77 

p-n 


w-Ct I 


u a i 

u-O f 

UQ t 

UC 1 


« © 1 ©O 9 

©D 1 

©r» i 

PO 1 

*-C 3 9 ©O 9 ©O 1 

►-© 



X2H | 

nez-. i 

« 7 t 1 

a7u | 

K 7 « 1 

Q £7 n 1 


» *5 1 1 X 2 " 9 

atz ft i 

at z n | 

otZ n 9 

•«z *» 1 •! 1 a z « 1 

atz ft 





< UJ | 

■•TUi t 

<UJ | 

«?U! 9 

•t 

1 © 1 — uj j 

0? UJ | 

< UJ | 

<tuj 9 

-t uj 1 < ui 1 <uj | 

fttUJ 

© 2 


P- Z 1 

P- ail 1 

p- 2 t 

*- >1 1 

U- Z 1 

© TT | 


9 l»- 21 

© IT 1 


© 7 1 

© 7 1 P IMP 71 

p- 7 

WO O 


WO o t 

to o 1 

wo O 1 

wo C2 1 

io © 1 

to c* 1 

40 

1 1 40 Dl 

> 1 t 
1 1 1 
* 1 1 

to O 9 

a i 

to 0 9 

WO O | 40 0 1 40 O 1 

1 I 1 

1 1 1 

1 1 I 

m» a 


i 


B8 



Table B2. Size Data (dn - cts/cm ) for 1977 Flights (continued) 
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Table B2. Size Data (dn - cts/cm ) for 1977 Flights (continued) 
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Table B3. Size Data (dn - cts/cm ) for 1978 Flights 
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Table B3. Size Data (dn - cts/cnr) for 1978 Flights (continued) 
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Table B 3 . Size Data (dn - cts/cm^) for 1978 Flights (continued) 
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Table B3. Size Data (dn - cts/cm ) for 1978 Flights (continued) 
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Table B3. Size Data (dn - cts/cm ) for 1978 Flights (continued) 
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Table B3. Size Data (dn - cts/cirr) for 1978 Flights (continued) 
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Table B3. Size Data (dn - cts/cnr) for 1978 Flights (continued) 
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Table B3. Size Data (dn - cts/cm ) for 1978 Flights (continued) 
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Table B3. Size Data (dn - cts/cnr) for 1978 Flights (continued) 
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Table B3. Size Data (dn - cts/cnT) for 1978 Flights (continued) 
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